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HAUPTAUFSÄTZE 


Beitrag zur T'heorie der langsamen Strömung zweier 
exzentrischer Kreiszylinder in der zähen Flüssigkeit. 
Von Wilhelm Miller in München. 


Im Anschluß an frühere Arbeiten des Verfassers wird die Strömung unter- 
sucht, die durch langsame Drehung zweier eine zühe Flüssigkeit beyrenzenden 
Kreiszylinder herrorgerufen wird und das System der Kräfte und Momente. 
die von der Flüssigkeit auf die Zylinder ausgeübt werden. Dabei wird Ge- 
legenheit genommen, einige bemerkenswerte Sonderfälle genauer zu dis- 
kutieren. 


1. Einleitung. 


In zwei größeren Arbeiten habe ich die langsame ebene Strömung untersucht, die durch 
die Drehung zweier Kreiszylinder erzeugt wird'). Insbesondere gelang es, neben der Strom- 
funktion die Zähigkeitskräfte in endlicher Form zu berechnen, die an jedem der beiden Zy- 
linder wirksam sind. Während die erste Arbeit von den Verhältnissen ausging, wie sie im 
allgemeinen bei einem im geschlossenen zylindrischen Lager rotierenden Zapfen vorausgesetzt 
werden, also vor allem die Wirkung auf den inneren rotierenden Zylinder untersuchte, bezog 
sich die zweite Arbeit im wesentlichen auf die Verallgemeinerung des Viskosimeter- 
Problems?), bei der der äußere Zylinder in Drehung versetzt und die Wirkung der Flüssig- 
keitskräfte auf den inneren (torsionsfähig aufgehängten) Zylinder bestimmt wird. In 
beiden Fällen war keine besondere Einschränkung über das Radienverhältnis und die Exzen- 
trizität der Zylinder angenomnren. Wie ich im besonderen in der ersten Arbeit gezeigt habe, 
lassen sich die Ergebnisse der Theorie ohne weiteres auch auf den Fall übertragen, daß beide 
Zylinder auseinander liegen oder der eine Zylinder in eine Ebene übergeht. In der vor- 
liegenden Arbeit möchte ich nun einige Ergänzungen geben, besonders mit Rücksicht auf die 
Annahme, daß beide Zylinder in gleichem oder entgegengesetztem Sinne in Umdrehung ver- 
setzt werden. Dabei wird es sich neben der Behandlung einiger bemerkenswerter Sonderfälle, 
für die übrigens auch die Strombilder berechnet wurden, in erster Linie darum handeln, das 
vollständige System der Zähigkeitskräfte aufzustellen und rechnerisch und graphisch zu er- 

1) Wilh. Müller: Ebene Spannungs- und Strömungsfelder mit zwei kreiszylindrischen Grenzen, Ing.-Arch. 
Bd. 13 (1942). 8. 37 bis 58. 

?2) Wilh. Müller: Über die Drehung von Kreiszylindern in der zähen Flüssigkeit und die Theorie der Kräfte 


am Rotations-Viskosimeter, Ann. Phys. Bd. 41 (1942), S. 335 bis 354. 
Beide Arbeiten, in denen auch die Literatur ziemlich vollständig angegeben ist, werde ich unter I und IT zitieren. 
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fassen, das in Verbindung mit beiden Zylindern auftritt, zunächst ohne Rücksicht auf die be- 
sonderen physikalischen Versuchsbedingungen. Bei der rechnerischen Behandlung, deren Er- 
gebnisse nur kurz gebracht werden können, verwende ich wieder die Methode der krumm- 
linigen bipolaren Koordinaten, die dem Problem am besten angepaßt ist und die Möglichkeit 
bietet, den ganzen Umkreis der Sonderfälle am vollständigsten zu entwickeln. 


2. Bipolare Koordinaten und Stromfunktion. 


Betrachtet man die beiden Kreise A, und K, als Elemente <=a, = eines Apollo- 
nischen Kreisbüschels £=konst. mit Beziehung auf die Pole P, und P, (= +ib), die zu- 
gleich die gemeinsamen Punkte des orthogonalen Kreisbüschels 7 = konst. sind, so läßt sich 
die komplexe Koordinate z=.x-+iy jedes Punktes der Ebene, bzw. des von beiden Kreisen 


begrenzten Feldes durch die bipolaren Koordinaten & und n in der Form darstellen 


bsinn bSin: 


2° CofE—cosy ""CofE— cosy 


Die Lage der Mittelpunkte und die Halbmesser der Kreise &= konst. und 7 = konst. sind dureh 


b 


(2) 

()=betgn, | 
bestimmt. Die Koordinate & hat im Pol P, (als Nullkreis) 
den Wert © und nimmt auf der positiven Seite der «- 
Achse positive Werte an, die abnehmen bis zum Werte 
Null, der dem unendlich großen Kreis £=0 (der x-Achse) 
entspricht. Die Koordinate 7 hat im höchsten Punkt A 
jedes &Kreises den Wert 7—=0, sie nimmt auf der posi- 
tiven Seite der y-Achse positive Werte, auf der anderen 
Seite negative Werte an und springt daher im tiefsten 
Punkt B jedes £-Kreises vom Wert +7 auf den Wert 
— 7 (Bild 1). 

Das hydrodynamische Problem läßt sich bei 
Beschränkung auf kleine Geschwindigkeiten oder kleine 
Reynoldssche Zahlen, bzw. Vernachlässigung der Träg- 
heitskräfte gegenüber den Zähigkeitskräften, zurück- 
führen auf die Lösung der biharmonischen Differential- 
gleichung für die Stromfunktion 


PB. Wenn wir bipolare Koordinaten einführen und Yin der Form 


Bild 1. Zur Einführung von b 
bipolaren Koordinaten und n. = Cof cos 


2 2 w 
ansetzen, so erhält man mit der Abkürzung 4*Y= an +5 


Für die Geschwindigkeitskomponenten in den Richtungen &£ und n erhalten wir 
1 


Der Rotor der (eschwindigkeit wird 
rtd=w=— 


| 
| 
{ | | 
| 
| 
XUM 
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Die Druckfunktion p genügt den Differentialgleichungen 


Von Konstanten abgesehen, bilden also der durch den Reibungskoeffizient « dividierte Druck 
und der Rotor den Real- und Imaginarteil einer Funktion wo=n+it der komplexen 


Variablen {=&+in, oder die Linien gleichen Druckes sind die orthogonalen Trajektorien 
zu den Linien konstanter Wirbelung. 


Was nun die Integration der Gl. (3) angeht, so läßt sich zeigen, daß die Stromfunktion 
unseres Problems in folgender Weise angesetzt werden kann: 


DEin2F)eosn + HESin:] . (8), 


wo die auftretenden Konstanten auf Grund der Randbedingungen und der Stetigkeitsforderung 
für den Druckverlauf ermittelt werden können. Wenn z. B. der innere Zylinder K, (£= a) 
mit der konstanten Geschwindigkeit », = — (v„);—.« rotiert, während der äußere = £ in Ruhe 


gs= 


bleibt, so ergeben sich durch eine Rechnung, die hier nicht wiedergegeben werden soll, bei 
der Festsetzung Y(ß)=0 für die Konstanten die Werte 


A, = [Einta— Pl: E,=2J, Coj’ Coj(a —P); 

B, F, =2J, Coj (a - G, =— J, (a Sin2ß; | 

2, Sin Sin’ (a 
NV 


D,—J, Cof (a+P): N 


wenn zur Abkürzung 


L,= (a — Sina — Sin Sin (a — | (9a) 


N = Ein (a — [la — P) (Sin? a + Sin? — 2 Sina Sin Sin (a — 
gesetzt wird. Damit ergibt sich 


(GofE 
N: = [(a — Sina — Sin(a— | 
(10). 


x P) (a — (Eof E— + Sin (E— PB) Eoj(E— a) cos] 
+ Sin Sin(a — Sin(E— P) - Sin Sin’(a— Sin 
— (a — Sina Cof — P) Sin 
Wenn dagegen der äußere Zylinder mit der Umfangsgeschwindigkeit — (vy):-3= v, rotiert 


und der innere ruht, so ergibt sich durch Vertauschung von « und £ und Ersetzung von ®, 
durch v, für die Stromfunktion der Ausdruck 


— cosn) 


N —= [Fin a Sin (a — PB) — (a — P) Sin X 

x —a) Eof (a — cosy) + Sin — a) Cof(£ — P) cos n] 
— Sina boj (a — PB) Sin (E— a) — Sina Sin’(a— a)SinE 
+ (a — Sin Coj a (a — P) Sin a) 


(il). 


Bei Drehung beider Zylinder hat man als Stromfunktion 


Bevor wir einige wichtige Sonderfälle untersuchen, wollen wir eine Übersicht über die Be- 
stimmung der Kräfte und Momente geben, die auf die einzelnen Zylinder wirken. 


3. Ergebnisse der Kräfteberechnung. 


Wir gehen am besten von dem Fall der Rotation des inneren Zylinders = a aus. Wie 


S 


ich näher begründet habe (vgl. I, S. 41), hat man für die Normal- und Tangentialspannung an 
einem Element des Zylinders die Ausdrücke 


—p; 


194? 
be- 
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Durch Einsetzung der Werte von ®,, und Benutzung von (8) ergibt sich 


(Eoja —cosn)[2J, Zin (a P)eosn+2J, Cof(la— . . . (1A). 


Aus der Form der Ausdrücke erkennt man ohne Schwierigkeit, daß bei der Integration über 
den Umfang des Kreises A, die y-Komponente der Kraft verschwindet. Es bleibt dann 


P,=-r,\peosgdg  r,\rsinpdg 

an 


Drückt man den Polarwinkel g durch , aus mit Hilfe der Formeln 


Cojacosn Sinady 
= 


so erhält man nach Auswertung der Integrale 
. ... . (6). 


In derselben Weise lassen sich die Kraft und das Moment berechnen, die auf den Zylinder K, 
übertragen werden. Dabei ist insbesondere auf den Vorzeichenwechsel zu achten, der dadurch 
entsteht, daß das Flüssigkeitsgebiet nieht mehr außerhalb, sondern innerhalb des Zylinders 
liegt. Man erhält dann 


Sina 
Daraus ergibt sich wegen 1,—r,<., .: 
Sin p 
Zinta- pP) 
H, Sin 
Sin (a — pP) 
Ein p 
Exzentrizität beider Kreise zum Radius des kleinen Kreises. Wir haben also 


Nun ist aber nach (2) der Faktor =: nichts anderes als das Verhältnis er der 


(18). 


M„+M,=-Pure= P,e,\ 
M„+M,+P,e=0 


Diese Gleichung und die Beziehung zwischen den Kräften stehen im Einklang mit dem Impuls- 
satz’), nach dem der Impulsfluß durch die Grenzflächen (Zylinder) bzw. sein Moment der Summe 
sämtlicher auf die Zylinder wirkenden Kräfte bzw. der Summe sämtlicher Kraftmomente 
gleich sein muß. Da der Impulsfluß verschwindet, so wird auch die resultierende Kraft bzw. 
das resultierende Moment gleich Null. 

Ganz entsprechend gestaltet sich die Rechnung für den Fall, daß der äußere Zylinder 
init der Geschwindigkeit », rotiert. Statt der Konstanten A,, B, usw. treten die durch Ver- 
tausch von a und 5 zu gewinnenden Konstanten A,, B, usw. auf. Man erhält dann die Werte 


=4nuH,; | 
„=4mur,(F, Sind + H, Coj 
— 4aur, (F,Sina+ a) | 
zwischen denen die Beziehungen bestehen: 


3) Wilh. Müller: Uber den Impulssatz der Hydrodynamik für bewegte Gefäßwände, Ann. Phys. 5. Folge 
Bd. 16 (1933), S. 480, 


Zina 

vos 

Coja 
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Setzt man nun die Werte der Konstanten ein, so entstehen schließlich folgende Formeln: 
Sin Sin(a—P 
13). 
(a — pP) (Sin’a+ Sin?) — 2 &in a Sin Sin(a— 
(14) (a — P) Sin’a Ctg (a — — (Sin’a — Sin’(a — P)) 
—P) (Sin’a + Sin? P) — 2&ina Sin Sin (a — 
- (Sim’a + Sin? p) 2 Sina Sin Zin(a- P) 
a Sin (a — pP) 
(a Sm’a+ -2L2ina P) 
(a Sin? (a - — — Sin’(a-- P)) 
(a — (Sin’a+ Sin? pP) — 2 Sin a Sin Sin(a— pP) 
Sin?’ [la — pP) Etg (a — p)— 1] 
— P) (Sin?a+ Sin?’ 2Sina Sin Sin (a — P) 


M.„=—tnur,v 


M„=4nrur,v 


Die Ausdrücke können noch vereinfacht werden dureh Einführung des Radienverhältnisses 
Sina Sin(a — 
und der Exzentrizität e= 
Sin p 


stätigenden Formeln: 


bzw. durch Verwendung der leicht zu be- 


ine; Sin (a — Lina; koja- 


(16). 

lurch 


ıders Cofla - = P= 


24 


2: Sina 
Wenn man gleichzeitig beide Fälle kombiniert, also voraussetzt, daß beide Zylinder mit den 
Geschwindigkeiten », =r,o, bzw. v,—=r,®, rotieren, so ergibt sich für die resultierenden 
Kräfte und Momente 


(1m). 


(a 


M,=4nur,r, 


Die zugeordneten Beiwerte kann man in verschiedener Weise abspalten. Setzt man etwa 


(18). 

P,=4rur,o,(,), =4rur, 
mme 1,0, (u), 
ıente 1 
bzw. so wird für »,=omw, bzw. w, m, 

N # 
2 
"erte (a —- P)(1+4) 2: Sina 


4 (1 |(a — P) Etg(a— P)- 1] 
(19), = (a— 2: Sina 
(m) (a— pP) 2: Sina 


(a— HB) 2: Sina’ (ep), = (a— pP) — 


Folge == 


Mech. 
(22). 
| 
0" 
= 
XYUM 


Z. angew. Math. Mech. 
Bd. 22 Nr. 4 Aug. 1942 
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Die Werte sind bereits für die Fäl- 
le, daß nur einer der Zylinder ro- 
tiert, eingehenddiskutiertund für 
verschiedene Radienverhältnisse 
inAbhängigkeit von der Exzentri- 
zität graphisch veranschaulicht 
worden (vgl. Tund II). Ergänzend 
haben wir hier (in den Bildern 2 
und 3) für ein bestimmtes Ra- 
dienverhältnis A=1,6 die Bei- 
werte (CH), (Cm)ı» (Em)s» 
(c',), für verschiedene Werte von 
o und 0’ in Abhängigkeit von © 
dargestellt. Es ergibt sich im 
besonderen, daß bei Drehung 
eines Zylinders die Momente an 
beidenZylindern entgegengesetzt 
gleich werden, wenn die Exzen- 
trizität gleich Null ist und daß 
bei wachsender Exzentrizität 
das Moment am ruhenden Zy- 
linder abnimmt, während das 
Moment am rotierenden Zylinder 
wächst. Für den theoretischen, 
physikalisch nicht verwertbaren 
Grenzfall der Berührung beider 
Zylinder (A—1=x=e) erhält man 
im allgemeinen unendlich große 
Kraft- und Momentenwerte. Nur 
an den ruhenden Zylindern ver- 
schwinden die Momente für e=x, 
d. h. (e/„), bzw. (c„), werden 
Null für 0 =0 bzw. o=0, oder 
beide Werte werden unbestimmt, | 
wenn 0’ bzw. o von positiven 
oder negativen Ausgangswerten 23 

gegen Null konvergieren. Bild 3. Kraftbeiwerte (cp)ı und (cp’)ı bei Rotation beider Zylinder. 


4. Bemerkenswerte Sonderfälle. 


Unter den vielen Spezialfällen, die möglich sind, wollen wir einige Fälle herausgreifen, 
die durch verhältnismäßig einfache Form der Stromfunktion gekennzeichnet sind. Wenn wir 
etwa die Kreise, also die Größen a und £ ungeändert lassen, so können wir z. B. das Ge- 
schwindigkeitsverhältnis so festsetzen, daß die Durchflußmenge dureh einen Querschnitt des 
Ringraumes zwischen den Zylindern einen vorgeschriebenen Wert annimmt, also z. B. ver- 
schwindet. Ferner lassen sich gewisse Sonderfälle durch bestimmte Annahmen über die Kraft- 
wirkung auf die Zylinder kennzeichnen. 


a) Eine besonders einfache Annahme ist das Verschwinden der Durchflußmenge 0 = (a) 
Y()=0. Da und P;(a)=V, so ergibt sich 
Mit Einführung des Nenners N (vgl. (9a)) ergeben sich nach einiger Rechnung folgende 
Ausdrücke: 


=», — B) Sina — Ein Sin (a Eof (a — p) — Ein (a — P)], 


(26). 


— B) Sin B -- Sin a Sin (a — B)]- — Cofla — Sin (a — P)] 


Daraus entnimmt man dann als Bedingung für die Drehgeschwindigkeiten 


(a — Sina — Sin (a—P) (a— — Sina 


Sin (a—P) — (a — B) Sinß "ırfe Sina—(a—P)] 


1 
| | 
06 
- 
XUM 
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oder 
(a— e Sina 
|e Sina — (a — P)] 


(a— pP) Sina — Ein Sin(a—P) 


Dir 
Zina Sina Sin (a— P) - (a— pP) 


= 


Wenn man ®, durch », ausdrückt und in 7, einsetzt, so erhält man für die resultierende 
Stromfunktion + Y, den Ausdruck 


- Pl - Sina Sin (a — P)] 
ev, (a pP) 
-(a P)’ Sina Sin’p Sin(£ -P)+ (a — Sin’a a) 
(a - B) Sina Sin? Sin (a — P) Ein(E a) — (a P) Sin’a Sin p Sin (a — PB) Zin(£— 

+(E a) Ein? (d — P) Sina Sin p P) Sin?’ (a — P) Sina Sin 

— (a — a) Sin’ a Sin(a - (a— P)Sin’PSin(a 
Dabei ist cos durch E=Lof (E — a) Coja+ Sin - a) Sina ersetzt, was mit Rücksicht 
auf den im Nenner auftretenden Faktor Coj £E— cosn gleichbedeutend ist mit der Abänderung 
der Konstanten. In dem Ausdruck Z der rechten Seite läßt sich aber der Nenner N heraus- 
heben. Es wird dann 
(a pP) EinpSin(E a) 


Zin(a pP) 


Dann entsteht für Y der relativ einfache Ausdruck 


-cosy) Ein(a— P) Ein! — (a - pH) Sind Sin (Su) (28) 
= Sin a Sin (a — PB) — (a — Sin 
Für die Konstruktion der Stromlinien kann man die reduzierte Stromfunktion 


a) Sin(a - P) — (a - pP) F(&) 


benutzen. Setzt man y'=e, so wird 


cos = „Fo. 


Jedem Werte von e entspricht eine Stromlinie. Läßt man e nach einer arithmetischen Reihe 
zu- oder abnehmen, so erhält man die Schar der Stromlinien derart, daß die Flüssigkeits- 
menge zwischen zwei benachbarten Linien gleich bleibt. Bei Übertragung auf rechtwinklige 
Koordinaten hat man 

besinn beSins 


- 29). 


Die Funktion y (bzw. Y) wird an den Kreisen X, und K,, d. h. für <=a und = Null. 
Die Stromlinie zerfällt also für e=0 in die beiden Kreise. Auf der Zentralen zwischen den 
Kreisen liegen nun zwei Punkte &, und £&,, in denen die Geschwindigkeit der Strömung Null 


wird. Sie ergeben sich aus 2; —0 für n=0V bzw. =r. Das ergibt die beiden Bestim- 


mungsgleichungen 


(a — Cofa — Sin(a— —a+ Ein!) - (a — P) Eina Sin 


(a — Ein B Cofa — Sin(a— PB) (&, a- (a - P) 


Wie die nähere Diskussion oder die Zeichnung ergibt, haben diese beiden Punkte ver- 
schiedenen Charakter. Während der im engsten Teil des Zwischenraumes gelegene Punkt 
S=£,,n—7 als Doppelpunkt einer ausgezeichneten Stromlinie angehört und eine Stelle be- 
zeichnet, an der die Strömung eine Umkehrung erfährt, ist der Punkt ©=&,, n=0 ein 
‚ isolierter Punkt oder der Mittelpunkt einer ihn umgebenden zirkulatorischen Strömung. Wir 
wollen den ersten Punkt als Umkehrpunkt (oder Staupunkt), den zweiten als Zirkulations- 


XUM 
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punkt bezeichnen. Die übrigen 
Stromlinien zerfallen in zwei 
Zweige, die durch die singuläre 
Kurve mit Doppelpunkt getrennt 
werden. 

In den Bildern 4 und 5 
haben wir die Stromlinienschar 
für zwei Fälle konstruiert, die 
durch folgende Werte als Haupt- 
daten gekennzeichnet sind: 


i. a=13%3, P= 05236: 
05286; 
1,= 18255; 

Umkehrpunkt 
—=0,9425, 
Zirkulationspunkt 
:,—=0,767, n=0(y,= 2,134 b). 


Singuläre Stromlinie entspricht 
c=c,=0,1%. 


2. a=1,669,  B= 1,0667: 
(9), = 1074b, 7, = 0,34 5; 
(Yo)s = 1,268 b, 0,7807 5. Bild 4. Stromlinienbild in der Umgebung zweier gegeneinander 
rotierender Kreiszylinder mit @=0, Fall I. 
2, —=13049, 7=0(y,=1,744b): 
ec, = 0,0383. 


Besonders einfach gestaltet sich 
die Berechnung der Druckfunktion. 
Man erhält 


Die Schar der Kurven gleichen 
Druckes, die für den zweiten der 
eben genannten Fälle in Bild 6 
wiedergegeben ist, stimmt übrigens 
im ganzen für alle Kreispaare 
überein, die demselben Werte von 
b zugeordnet werden können und 
für de ist. 


el). 
(a — P) 
Zina pP) -(a— 


Für die Kräfte und Momente 
erhalten wir die Ausdrücke 


p) 
M,=-4rurio, 
(a-p) 
M,=4aur}o, x 


e&oj 


FR. rg _ Bild5. Stromlinien in der Umgebung zweier rotierender Zylinder 
Sina — (a— P)] mit Q=0, Fall I. 
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Zwei Unterfälle sind noch be- 
merkenswert: 

a) Wenn 3=0 wird, so geht 
der Kreis K, in eine im Unendlichen 
als geschlossen zu betrachtenden ge- 
rade Linie über, die entsprechend der 
Gl. (27) mit der Geschwindigkeit 


a 


(33) 
in sich selbst verschoben wird. Die 
Stromfunktion nimmt die einfache 
Form an 


(34). 


Für den Doppelpunkt (Umkehrpunkt 
der Strömung) hat man mit „= 


Sinö+i-a=0. . . (5). 
Der Zirkulationspunkt rückt hier ins 
Unendliche. In dem gezeichneten 
Beispiel (Bild 7) ist a=1,647, 

Bild 6. Kurven gleichen Druckes bei zwei gegeneinander 
rotierenden Zylindern (Q=0). = 0,1825, 0,3724 b. 
Für die singuläre Stromlinie ist 


1 
0,32, 
v, 


P) Wenn negativ 
wird, so liegt der Kreis K, 
auf der negativen Seite 
der Achse &=0 und ganz 
außerhalb des Kreises 
Die entsprechenden For- 
meln werden besonders 
übersichtlich, wenn «a, 
also beide Kreise einander 
gleich sind. 


Die Bedingung (27) 
liefert in diesem Falle 
v,=v,=r. Beide Kreise 
müssen also im gleichen 
Sinne mit gleicher Ge- 
schwindigkeit sich drehen, 


it Q=0 ist (s. Bi auf- 
Bild 7. Strömung um einen rotierenden Kreiszylinder in der Nähe einer tangential damit Q x ist (s Bild8a f 
bewegten ebenen Wand mit Q=0. S. 188). Die Stromfunktion 


geht über in 
2bv Coja-ESin?— adinakoj 


GofE—cosy 
Die Stromlinien, die bereits früher gegeben wurden (vgl. I, S. 48), erinnern an die sog. 


Cassinischen Kurven. Der Doppelpunkt rückt in die Mitte der Zentralen beider Kreise. 
Die weitere Diskussion ist in der angeführten Arbeit I nachzulesen. 


b) Der Fall verschwindender Kraftwirkung. Die Kraft, die auf einen der 
Zylinder wirkt, verschwindet, wenn 


oder die Momente der Umfangsgeschwindigkeiten in bezug auf die Mittelpunkte entgegen- 
gesetzt gleich werden. Wir können also setzen 


A; Nr.4 Aug. 1992 Müller, Beitrag zur Theorie der langsamen Strömung zweier Kreiszylinder 187 


Durch Einsetzung in die zusammengesetzte Stromfunktion Y=%Y,+Y, erhält man einen 
Ausdruck von der Form 


Sina: _ N ElEofE — cosn) — Zinsi+B, 
+ (a — Sin’a Sin - P)EojE+(a — Sina Sin? Coja Cof (a Cof & 
Sind (a - [Zin (a + P)Coj(E - a)cosy 
+ Sin (£ — a) Eof cos n] 
gesetzt ist.. Den Ausdruck B kann man umgestalten mit Hilfe der Beziehungen 
Sin — a) + Sin (E — a) Eof ($ H)— Zin(a+P) 
—=2 Lin Sinla+P). 
Dann läßt sich die Konstante N abspalten, und es wird 


Sin (E — pP) 


Schließlich erhalten wir nach einiger Rechnung die Stromfunktion 


Sin (a — P) (Eof E— cos n) 


v,r, Gig (a = P) 


Die Durchflußmenge wird 
9= Pla) - P&tgla- M- 1] 
Für die auf die Zylinder wirkenden Momente ergeben sich die Werte 
M=-M, 2 2 2 22. (40). 


a) Erwähnenswert ist zunächst der Sonderfall 5=0, für den v,—=0 wird. Wir er- 
halten also den Fall der Rotation eines Zylinders in der Nähe einer ruhenden ebenen Wand 
mit der weiteren Bedingung, daß auch im Unendlichen Ruhe herrschen soll. Dann wird 


Sin a — a)cos 


Sin a (Cof cos n) 


Das Bild 9 auf Seite 188 zeigt den Verlauf der Stromlinien. Die ausgezeichnete Strom- 
linie zerfällt in die Gerade £&=0 und eine geschlossene Kurve, die £=0 berührt. Der Doppel- 
punkt wird also hier zum Selbstberührungspunkt. Für das widerstehende Moment) findet man 


Wenn 5 unendlich groß wird, also die Wand ins Unendliche rückt, so erhalten wir als 
Reibungsmoment eines rotierenden Zylinders in der unendlich ausgedehnten Flüssigkeit 


P) Wenn = -a wird, so erhält man &,=—m,. Die gleichen Kreise rotieren ent- 


gegengesetzt (Bild 10). Dann wird die Stromfunktion und das widerstehende Moment 


% Vgl. auch Wilh. Müller: Einführung in die Theorie der zähen Flüssigkeiten, Leipzig 1932, S. 220. 
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ICH 


Bild 11. 


Rotierende konzentrische Zylinder. 


Bild 8. Kräfte und Momente bei zwei gleich 
sinnig rotierenden gleichen Zylindern. ° 


Bild 10. Gegensinnie rotierende gleiche 
Bild Ss. Bild In. Zylinder. 


Bild 9. Strömung um einen rotierenden Kreiszylinder in der Nähe einer ruhenden ebenen Wand. 


Zins [Coj’a 


« 
a+ ZinaCoja 


(45). 


c) Konzentrische Zylinder. Wenn die Zylinder konzentrisch sind, so lassen sich 
die entsprechenden Ausdrücke z. B. für das Moment durch Grenzübergang bestimmen. Wir 
haben zunächst 


(€ in al Coj = (A? Zin (a - P) = 


1 77 
© 
| 
| 
| Ach 
| 
| 
= 
XUM 


4. angew. Math. Mech. Müller 
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2+1 —1 
(a— p)Etala - -1= Ar Tg 


(—1). 
Damit erhält man für das auf den Zylinder wirkende Moment 


4 
M=-M,=-Anur,r,(o, Un Er 


5. Bemerkungen zur Lösung der ungekürzten Differentialgleichung. 


Wenn man die Trägheitsglieder in der hydrodynamischen Differentialgleichung berück- 
sichtigt, so erhält man als Gleichung für die Stromfunktion Y der stationären Strömung zu- 
nächst in rechtwinkligen Koordinaten 


dw 


Bei Einführung von bipolaren Koordinaten geht die Gleichung über in 


oder ausführlicher 


— Sın . 2 4 
0 07 05 


Um eine weitere Annäherung an die exakte Lösung dieser Gleichung zu erhalten, kann man 
z. B. für den verallgemeinerten Zapfenfall O< <a, v,=0), also bei Rotation des inneren 
Zylinders die Abweichung der Stromlinien von den Kreisen &—= konst. als klein vernachlässigen 


und «daher = 5 = setzen. Dann lautet die vereinfachte Differentialgleichung 


vn 


In (48). 


Die Lösung dieser Gleichung dürfte sich für große Reibung am zweckmäßigsten in der Form 


einer nach Potenzen der Reynoldsschen Zahl R,= gi we fortschreitenden Reihe 


ansetzen lassen, in der Y, die erste für große Zähigkeit gültige Näherung darstellt, wie sie 
oben benutzt wurde. Die Bestimmung der Funktionen %,, 7,, ..., die mehrere Schwierig- 
keiten bietet, ist bisher nicht in Angriff genommen worden. Wenn sie gelingt, etwa für den 
Fall kleiner Exzentrizität und kleiner Radiendifferenz, dürfte es auch möglich sein, über die 
zweite, in die Zentrale der beiden Kreise fallende Kraftkomponente, die auch durch Versuche 
festgestellt wurde, nähere Angaben zu machen und damit z. B. die Theorie der Kräfte am 
Zapfen weiter zu vervollständigen. 364 
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Über die Luftkräfte, die auf einen 
harmonisch schwingenden zweidimensionalen 
Flügel bei Überschallgeschwindigkeit wirken. 

Von S. vr. Borbely in Kolozsvär (Ungarn). 
Bericht des Flugtechnischen Instituts der Technischen Hochschule Berlin. 


Die Arbeit enthält einen Beitrag zur Theorie der Strömungen mit Über- 
schallgeschwindigkeit. Es wird das Problem des harmonisch schwingenden 
ebenen Flügelprofils behandelt, indem das Strömungsfeld aus gewissen 
Elementarlösungen der linearisierten Differentialgleichung. die mit Hilfe 
ron Hakenintegralen gewonnen sind. aufgebaut wird. 


I. Übersicht. 


Vorliegende Arbeit behandelt die Strömung um einen schwingenden Tragflügel, der 
mit Überschallgeschwindigkeit angeströmt wird. Es werden für den Fall des ebenen Problems 
(Flügel von unendlich großer Spannweite) unter Voraussetzung von Reibungsfreiheit die Luft- 
kräfte ermittelt, die bei harmonischer Dreh- und Biegeschwingung des Flügels entstehen'). 

Die mitgeteilte Theorie umfaßt alle Fälle, bei denen die vom Flügel herrührendeu 
Zusatzgeschwindigkeiten «, » an allen Stellen klein gegen die Anströmgeschwindigkeit U, sind. 
Praktisch trifft das zu, wenn folgende Bedingungen erfüllt sind: 


1. Der Anstellwinkel des Flügels sei klein, die Schwingungsamplitude sei klein, und 
die Schwingung erfolge langsam (diese Bedingungen treten in gleicher Weise bei 
der Theorie der Flügelschwingungen in inkompressiblem Medium auf; siehe z. B. 
Kassner-Fingado, Das ebene Problem der Flügelschwingung, Luftf.-Forschg. 
Bd. 13 (1936), S. 374 bis 387. 

2. Das Flügelprofil sei scharfkantig und schlank, so daß die Neigung der Tangenten an 
die Profilkontur überall klein ist. 


Es läßt sich dann zeigen, daß die Strömung, die durch einen mit Überschallgeschwindigkeit 
angeströmten schwingenden Flügel erzeugt wird, in zwei Anteile zerlegt werden kann’): 


I. einen stationären Anteil; dieser stellt die Strömung um das gegebene, mit dem 
mittleren Anstellwinkel der Schwingung gegen die gegebene Anströmgeschwindigkeit 
fest angestellte Profil dar. Der stationäre Anteil ist also durch die Form, d.h. 
Wölbung und Dicke, und den Anstellwinkel des Profils bedingt. 


. einen instationären, harmonisch schwingenden Anteil; dieser stellt die Strömung um 
eine ebene, symmetrisch gegen die Strömung schwingende Platte von derselben Tiefe 
und derselben Schwingungsfrequenz und -amplitude wie der gegebene Flügel dar. 
Dieser Anteil ist also von Form und Anstellwinkel des Profils unabhängig. 


Die auf den Flügel wirkende Luftkraft erhält man durch Superposition der Kräfte aus 
beiden Strömungsanteilen. Der erste, stationäre, Strömungsanteil ist durch die Arbeiten von 
Prandtl, Ackeret, Busemann bereits bekannt. 

Den Hauptinhalt der vorliegenden Arbeit bildet die Untersuchung des zweiten, in- 
stationären, Strömungsanteils. Die Luftkräfte, die sich bei dieser Strömung ergeben, findet 
man in Gl. (31) angegeben. Die dort auftretenden Glieder g,, r; sind Funktionen der 
Machschen Zahl und einer reduzierten Frequenz, die sich nicht elementar angeben lassen. 
Die einzelnen Glieder der Luftkraft lassen sich als Massen-, Auftriebskraft usw. deuten, wie 
es in ähnlicher Weise auch bei Schwingungen in inkompressibler Flüssigkeit möglich ist. 
Um zu zeigen, wie sich die Kräfte bei Übergang von inkompressibler Flüssigkeit auf Über- 
schallströmung ändern, sind in Tafel I die Luftkraftanteile bei Überschallströmung und bei 
inkompressibler Flüssigkeit (nach den Angaben von Kassner-Fingado) gegenübergestellt. 

Die mathematische Behandlung des Problems geht aus von der Linearisierung der 
Differentialgleichung. Mit Hilfe von Hakenintegralen bekommt man das Strömungsfeld eines 
Druckpunktes, der sich stationär oder instationär (harmonisch schwingend) bewegen kann. 
Aus derartigen Druckpunkten läßt sich das Strömungsfeld um ein Profil aufbauen. 

Zum Schluß wird gezeigt, wie sich auch die stationäre Strömung mit Hilfe der Haken- 
integrale behandeln läßt. 

1) Nach Abschluß der Arbeit wurde mir bekannt, daß folgende Arbeit sich mit dem gleichen Thema befaßt: 
Possio, L’azione aerodinamieco sul profilo oseillante alle veloeitä ultrasonore. Acta Pontif. Acad. Sei. Bd. 1 (1937). 


Nr. 11, S. 9. 
?) Diese Zerlegung ist ganz analog der bei inkompressibler Flüssigkeit möglichen. 
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Vorliegende Arbeit ist auf Anregung des damaligen Leiters des Flugtechnischen Instituts 
Prof. Dr.-Ing. Herbert Wagner entstanden. Für seine Ratschläge sei ihm auch an dieser 
Stelle gedankt. Besonderen Dank schulde ich ferner Dr.-Ing. H. Simon, der sich der Mühe 
unterzogen hat, nicht nur die ganze Arbeit kritisch durchzusehen, sondern sie auch für die 
Veröffentlichung zu überarbeiten. 


II. Bezeichnungen. 


em Ortskoordinaten 
see Zeit 
em Koordinaten der Drehachse des Flügels (Bild 1) 
Anstellwinkel des Flügels, d. h. der Profilsehne (Bild 1) 
em Amplitude der Schlagschwingung 
Amplitude der Drehschwingung 
Schlag- und Drehfrequenz 
Flügeltiefe 
die von der Profilsehne aus gemessene Ordinate der Profilkontur (Bild 1) 
em sec ' Geschwindigkeitskomponenten der Flüssigkeit in «- bzw. y-Richtung 
em sec"' Komponenten der Störgeschwindigkeit in «- bzw. y-Richtung 
em see" ' Anströmgeschwindigkeit der Flüssigkeit 
em? s Geschwindigkeitspotential der Störgeschwindigkeiten 
em see” ' Schallgeschwindigkeit 
kg sec? em”* Dichte 
kg em”? Druck 
Verhältnis der spezifischen Wärmen bei konstantem Druck und konstantem 
Volumen 
Auftriebskraft 
Moment der Luftkräfte um die Flügelvorderkante, im mathematisch positiven 
Sinne 


sin Machscher Winkel 


ot 


reduzierte Frequenz 
cos? 


Querstriche («, usw.) bedeuten im mathematischen Teil dimensionslose Variable (Koordinaten). 


Indizes: jı= am Flügel 
o, u oben, unten 
0 bezieht sich auf die ungestörte Strömung vor dem Flügel. 


III. Voraussetzungen und Annahmen. 


Wir betrachten eine reibungsfreie unendlich ausgedehnte zweidimensionale Flüssigkeit 
(ebenes Problem), die vom Unendlichen her mit gleichmäßiger Geschwindigkeit gegen einen 
Flügel anströmt; vom Einfluß der Schwere auf die Strömung wollen wir absehen. Die Ge- 
schwindigkeit der Anströmung sei größer als die Fortpflanzungsgeschwindigkeit des Schalls 
in der betreffenden Flüssigkeit. 

Um die Behandlung der Strömungsvorgänge zu vereinfachen, betrachten wir den 
mathematischen Grenzfall, daß das gleichmäßige Geschwindigkeitsfeld der Anströmung durch 
den Flügel an allen Stellen nur sehr (unendlich) wenig geändert wird, d. h. daß die durch 
den Flügel erzeugten Zusatzgeschwindigkeiten «, » sehr (unendlich) klein gegen die An- 
strömgeschwindigkeit U, sind. Dazu ist erforderlich, daß 


I. der Anstellwinkel des Flügels sehr (unendlich) klein ist und daß die Bewegungs- 
geschwindigkeit des Flügels in Richtung senkrecht zur Anströmgeschwindigkeit sehr 
(unendlich) klein gegen die Anströmgeschwindigkeit ist°). 


2. die Neigung der Profiltangenten gegen die Anströmrichtung überall sehr (unendlich) 
klein ist, d. h. also, daß das Profil scharfkantig und schlank ist. 


In einer Überschallströmung ist die Zusammendrückung des strömenden Mediums für 
die Strömungsvorgänge wesentlich. Diese Zusammendrückung ist im Gebiet der Strömung 
im allgemeinen ein stetiger Vorgang, d. h. Druck p und Dichte o ändern sich von Stelle zu 
Stelle allmählich, und zwar adiabatisch nach der Gleichung p=a:o*. Der Vorgang verläuft 


3) Es sei bemerkt, daß diese Voraussetzungen auch bei der Tragflügeltheorie in inkompressiblem Medium 
üblich sind. 
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ohne Entropievermehrung und ist daher umkehrbar. Es kann jedoch gewisse Unstetigkeits- 
stellen geben, an denen „Verdiehtungsstöße* auftreten, d. h. an denen Druck und Dichte 
piötzlich von einem Wert „auf einen Wert , übergehen. Da diese Verdiehtungsstöße nicht 
adiabatisch, sondern unter Entropievermehrung verlaufen, sind sie nieht umkehrbar. Es gilt 
für sie die Gleichung 


intwickelt man nun das Dichteverhältnis o,/o, für Stoßvorgänge und adiabatische Vorgänge 
nach Potenzen von (p, — P,)/P, = Ip/Pp,., So erhält man in beiden Fällen: 


Po Po 


d. h. beide Verdiehtungsvorgänge stimmen bis zu Gliedern zweiter Ordnung einschließlich 
überein. Da in dem von uns betrachteten Grenzfall sehr kleiner Zusatzgeschwindigkeiten 
die Glieder höherer Ordnung gegen die linearen Glieder verschwinden, so besteht kein Unter- 
schied zwischen adiabatischer und Stoßverdiehtung; wir können also überall adiabatische 
Zustandsänderungen annehmen. 


Aus der eben angestellten Betrachtung folgt ferner, daß wir die Strömung im ganzen 
Gebiet als rotationsfrei-ansehen dürfen. Bewegt sich nämlich mit dem strömenden Medium 
eine beliebige geschlossene Kurve, so gilt für die zeitliche Änderung der Zirkulation /' längs 


dieser Kurve 
d t oO 


wobei das Integral längs der Kurve zu nehmen ist. Treten nur adiabatische Zustands- 
änderungen in der Strömung auf, so verschwindet dieses Integral; sind auch Stoßvorgänge 
vorhanden, so verschwindet es nicht mehr, aber seine Größenordnung ist nach dem Voran- 
gehenden (1p/p,)‘. Die Zirkulationsänderungen verschwinden also in dem von uns betrachteten 
Grenzfall, und da die Parallelströmung vor dem Flügel rotationsfrei ist, muß die Bewegung 
im ganzen Strömungsgebiet wirbelfrei bleiben. 


Schließlich sei noch vorausgesetzt, daß der Flügel eine harmonisch schwingende Dreh- 
und Schlagbewegung senkrecht zur Anströmrichtung ausführe, und zwar schon seit sehr 
(unendlich) langer Zeit, so daß hinter dem Flügel eine in «-Richtung unbegrenzte periodische 
Strömung vorhanden sei. 


IV. Die linearisierte Differentialgleichung. 


Nach unseren Voraussetzungen lassen sich die Geschwindigkeiten U, V, der Druck p 
umd die Dichte o in folgender Form darstellen: 


7) 0, (2%, y; 7) | 
dl), 


wobei u/U,, r/U,, 0,/0,, p,/p. Kleine Größen sind. 


Für den Zusammenhang zwischen Druck und Dichte p= «a: 0“ läßt sich dann schreiben: 


p=ao0,taxo 


d —)axo,” ' die zur ungestörten Strömung’, gehörige Schall- 


Dabei bedeutet ce — 


geschwindigkeit. 


Führen wir diese linearisierten Ausdrücke in die Kontinuitätsgleichung 


do, 


und die beiden Eulerschen Gleichungen 


-Scheel, Handbueh der Physik, Bd. VH, Berlin 1927, S. 325. 


4, Siehe z. B. Ackeret, Gasdynamik in Geiger 
5) Siehe Ackeret, a. a. O., S. 294. 


»—1\o, 2 0 9ı 


Kin v. Borb&ly, Über die Luftkräfte bei Überschallgeschwindigkeit 


1op 
dx 


ein und vernachlässigen Glieder zweiter und höherer Ordnung, so gehen die Differential- 
gleichungen über in 


(4) 


ou do 
— 
0,08 


Eliminiert man aus diesen Gleichungen o,, so erhält man 
woraus schließlich mit Hilfe der Bedingung der Wirbelfreiheit 
du 
für die y-Komponente der Differentialgleichung folgt: 


+21 


V. Randbedingungen; Zerlegung des Problems. 


Wenn wir die Randbedingungen, denen die Lösung der gewonnenen Differential- 
gleichung (6) genügen muß, mathematisch formulieren, so ergibt sich in unserem Falle: 


1. Im Bereich vor dem Flügel muß, da in einer Überschallströmung sich Störungen nur 
stromabwärts auswirken können, die ungestörte Parallelströmung vorhanden sein, 
d. h. für 

2. Im Bereich des Flügels muß die Strömung am Flügel anliegen, d.h. in der Normal- 
richtung zum Flügel muß die Geschwindigkeitskomponente der Strömung mit der des 
Flügels übereinstimmen. Wegen der Voraussetzung der (unendlich) kleinen Neigung 
der Flügelkontur ist die Normalkomponente mit der Vertikalkomponente V identisch. 
Wenn die Lage des Flügels in der Form y=g(#,r) (im Bereich 0....x...t) gegeben 
ist, so muß also am Ort des Flügels, d. i., bei unserer Beschränkung auf die linearen 
Glieder, für y=0 (und O...x...f) gelten 

dg_ 09, 

3. Im Bereich hinter dem Flügel (#« >t, y=0) muß die Strömung entsprechend der 
Schwingungsfrequenz des Flügels periodisch (zeitlich) wechseln. 


Wir wollen nun zeigen, daß das Strömungsproblem in zwei Teil- 
probleme, ein stationäres und ein instationäres, zerlegt werden kann. 


Wir betrachten die Schwingung des Flügels. Sie sei als Summe 
(2). N einer Schlagschwingung und einer Drehschwingung um eine Achse mit 
| der Abszisse x, (siehe Bild 1) gegeben. Die Ordinate eines beliebigen 
rall- Flügelpunktes mit der Abszisse « hat, wenn wir nur die Glieder 1. Ord- 
nung berücksichtigen, den Wert: 


- 


Io,u(&) = Ya — (r — ul) 


Der Ausschlag y. und der Drehwinkel e sind dabei harmonische (gegeneinander phasen- 
verschobene) Funktionen der Zeit mit der Schwingungsfrequenz o: 


Y=Yıt Ye‘) 
| 


(9). 


%) Wir verwenden die bekannte komplexe Schreibweise für die harmonische Schwingung; will man auf die 
reelle Schreibweise übergehen, so hat man entweder den Real- oder den Imaginärteil für sich zu nehmen. 


13 
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Führen wir diese Ausdrücke in die Randbedingung 2 ein, so lautet sie: 
für 0...x...t, 


Hierin stellt der erste Teil 


Ale 


dr 


die Geschwindigkeitsverteilung an einer um die Nullage schwingenden unendlich dünnen 
ebenen Platte dar, der zweite 


die Geschwindigkeitsverteilung an dem gegebenen, aber in der Mittellage fest angestellten 
Profil. Es liegt daher nahe, an eine Zerlegung des Problems in die beiden Teilprobleme der 
schwingenden Platte und des stationären Profils zu denken. Da die Differentialgleichung für 
v linear ist und somit das Superpositionsprinzip gilt, so steht einer Zerlegung nichts im Wege, 
falls auch die übrigen Randbedingungen dabei erfüllt werden können. Nun ist in beiden 
Fällen im Bereich vor dem Flügel v»—=0; hinter dem Flügel ist bei stationärer Strömung 
ebenfalls »—=0, während bei der schwingenden Platte dort die Strömung periodisch sein 
muß. Man sieht also, daß bei Überlagerung der beiden Strömungen wirklich alle Rand- 
bedingungen der Gesamtströmung erfüllt werden. Die Zerlegung ist also zulässig. 


Wir beschäftigen uns zunächst und hauptsächlich mit dem Problem der dünnen, nicht 
angestellten, schwingenden Platte. 


VI. Die Elementarlösung und der Lösungsansatz. 
Eine Partikulärlösung der Differentialgleichung (6) ist die Funktion 


eirt+iuxtivy (g,u reell), 


wenn » die aus der charakteristischen Gleichung folgende Bedingung erfüllt: 


= no u . 


Da die Differentialgleichung linear ist, darf man die Lösung mit einem beliebigen, von ., y, t 
unabhängigen Faktor multiplizieren und verschiedene Lösungen überlagern. So kann man 
aus der Partikulärlösung folgende Elementarlösung bilden 


ve. 
“+7 
o 


Die Hakenintegrale, die in dieser Lösung auftreten, haben die Besonderheit, sich an gewissen 
Stellen diskontinuierlich zu ändern. So ist z. B. längs der «- Achse’): 


für x,: v,(2,0;)=0 


(—Xı) 


d. h. auf der «- Achse herrscht links von der Stelle «, die ungestörte Parallelströmung, während 
rechts von &, eine zeitlich harmonische Störgeschwindigkeit » 
vorhanden ist. Untersucht man den ganzen Bereich der .x y-Ebene, 
so findet man, daß der Bereich, in dem die Geschwindigkeit » 
von Null verschieden ist, einen Keil mit der Spitze in .«, darstellt 
(Bild 2). Die Elementarlösung stellt also eine gewisse Störung der 
Parallelströmung dar, die vom Punkt x, ausgeht. 

Um die Strömung um die schwingende Platte darzustellen, 
benutzen wir eine Verteilung derartiger Elementarstörungen über 
die Flügeltiefe, und zwar machen wir folgenden Ansatz: Bild >. 


?) Siehe darüber auch Abschnitt VII. 


€ 


ıaAA 
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eiuxtiry 
du+l(, \ du 


eier eiuxtivy 


— — 


t 
\ u+ u, 


>> xı= 
— 


t 
dud., - dud., 
Us 


wobei 
6 


> 
. (12a) 
Die Stärke der Störung an der Vorderkante, ©, bzw. C,, sowie die Verteilungsfunktion f («,) 
sind dabei vorläufig noch unbestimmt. Wie leicht ersichtlich °), liefern für y=0 die Haken- 
integrale 


Rs), 


w 


und man erkennt, daß auf der «-Achse vor dem Flügel »—0, hinter dem Flügel » zeitlich 
(und örtlich) periodisch ist. Die Bandbedingungen 1 und 3 des $ V ur daher erfüllt. Es 
fehlt nun noch die Randbedingung 2, die im Bereich des Flügels 0. .t gilt. Schreibt 
man sie für die in Gl. (10a) gegebene Schwingungsform des Flügels auf und berücksichtigt, 
daß sich für y„=0 die Geschwindigkeit » der G]. (12) wegen (13) vereinfacht, so erhält man 


| 


x 

+ \ 


Diese Bedingung läßt sich mit dem Ansatz 


fa@)= a + 


erfüllen; sie geht damit über in: 


U 
x.) 
20 
e-ins+ P o-imx a 


+! 
iu, (iu im 


+C, |, — e-ilex 


iu, im, iu, 


d. h. die Konstanten C,, (,, a, ö müssen die Bedingungen erfüllen: 


1 1 r —ior 


Be-ion, 
2 


iu, 
Ihre Größe bestimmt sich daraus zu 
8) Siehe z. B. Rothe-Ollendorff-Pohlhausen: Funktionentheorie. Berlin (1931); Breisig: Theoretische 


Telegraphie, 2. Auflage. Braunschweig 1924. 
13° 
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0 


a= — 
31 


Y | Ee-ie 


B=—a7r Beric" 


= 


io 


- c) 


Y- ( 


e E 
+57, 


wobei zum Zwecke einer späteren Umformung hier schon die Gleichung vermerkt sei: 


„te 


(16). 


. r U, ior dg, 


Nachdem die Konstanten a, f, C,, C, bestimmt sind, ist die Geschwindigkeit v (x, y; 7) 
(Gl. 12) nunmehr bekannt. Wir wollen sie jedoch in etwas anderer Form schreiben, da dies 
zur praktischen Auswertung der Hakenintegrale zweckmäßiger ist, Wir führen eine neue 
Integrationsveränderliche » ein durch die Gleichung 


. (17a) 


und setzen 
. (17b\. 


>» 

— Ha —1, +1 


Führen wir weiter im äußeren Integral des Ausdruckes für » die dimensionslosen Größen 


ein, so erhalten wir schließlich für die Geschwindigkeit die Gleichung 


U 


—ior 


= 


VII. Zurückführung auf Besselfunktionen. 


In dem Hakenintegral 
eikw+igVw: —ı 
y(z,9;+D= \ o+1 
(— 1,0) 


ist der Integrand eine mehrdeutige Funktion von » °); die Riemannsche Fläche des Inte- 
granden ist zweiblättrig, im Unendlichen schlicht und besitzt die Verzweigungspunkte (+ 1,0). 


%, Es sei bemerkt, daß wir uns im folgenden auf positive Werte y, d. h. die obere Halbebene der Strömung, 
beschränken. 


196 | 
| 
— 
oc 
Ts 
Dann können wir schreiben 
+7 
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sro v. Borb&ly, Über die Luftkräfte bei Überschallgeschwindigkeit 


Wir führen den Integrationsweg H, der sich über die 
ganze reelle Achse erstreckt, so, daß die Verzw eigungs- 
punkte nach Bild 3 umgangen werden. Wählen wir 
den Zweig der Wurzel, der negativ wird, wenn @ 
reell und >1 ist, so geht der Exponent im Inte- 
granden mit wachsendem Absolutwert von 
asymptotisch gegen den Wert 


Rising -i-.cosp); 


der Realteil dieses Ausdrucks ist im Falle «— y„>0 
in der oberen »-Halbebene, d.h. für >0, nega- 
tiv, im Falle & -y<0 in der unteren Halbebene. 
Daraus folgt, daß das Integral 
| +iy Vor — 1 
über den Halbkreis K, der oberen bzw. K, der unteren Bild 3. 
Halbebene (für >0 bzw. 2 y< 0, siehe Bild 3) 
im Grenzfalle R>» gegen den Wert Null geht. Wir können das Hakenintegral daher durch 
Zufügung des oberen bzw. unteren Halbkreisintegrals zum geschlossenen Umlaufintegral er- 
gänzen, ohne seinen Wert zu ändern. 
Da im Fall « —y< 0 im Innern des Integrationsbereiches keine Singularität liegt, ist 


nach dem Cauchyschen Satz der Wert des Umlauf- und damit auch des Hakenintegrals 
gleich Null. 


1. für y<0 ya; [ =1im ( [+ 6. 


Im andern Falle dagegen liegen im Innern des Integrationsbereiches Singularitäten; da 
diese nur auf dem von (— 1,0) bis (+ 1,0) laufenden Verzweigungsschnitt liegen, dürfen wir 
den Integrationsweg auf den Umlauf S um den Verzweigungsschnitt zusammenziehen (Bild 3), 
ohne daß dadurch der Wert des Integrals geändert wird. 


2. für x -y>0 +f)= 


Wir zerlegen nun den Integrationsweg in zwei Teile; den einen Teil bildet ein Kreis 
vom Radius e um den Punkt (-- 1,0), den anderen Teil bildet der Rest des Umlaufweges. 


v(%,9,+ 1)= lim —— ie el!dp + deo|- 


Der erste Teil stellt das 2i-fache Residuum des Integranden am Punkte (-- 1,0) dar und 
hat den Wert 2rie-i*x; der zweite Teil ist der Cauchysche Hauptwert des Integrals. 


+l)=2nie-i!+ 
—1,+1,—1 
Der Hauptwert läßt sich durch die Substitution ® =cos t überführen in: 


eifcost + Teint 
—1,+1,—1 


ei iy Vo? —1 


sintdt. 


eix cost) eiYy(l+ cost) 
\ 


t 
1+cosi ei 1+cost 


ist, so erhält man: 


«w-£bene 
K, 

le 

(Ka) 
(—1,0) (— 1,0) 
Da nun 

0). 

ug, 

YııM 
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x +7 

= eiscost +Ysint .sintdtd&+ sint-dt. 


Das innere Integral des ersten Gliedes hat den Wert '"): 


iS cos Ysi 5 /=2 =s 
\ eiseost +Y y)=2ny 


wobei J,, /, Besselsche Funktionen erster Art bedeuten; damit ist: 


eizu -iyVur—1 —1 — 
-do= AUCH Mas 


+1 
—1,+1,—1 


(207. 


+ 7 
eiyieost—isint 


e-i@-n sint- dt 


2 


Das zweite Integral kann durch die Substitution eos !— = | 


z-Lbene oder z=cost—i-sint in ein Integral über einen (offenen) 
Kreis X der z-Ebene (Bild 4) umgewandelt werden: 
+7 
(10) +70) ei (cost isind) eidt.(z 1) 
T 


(K) 


Da der Integrand Residuen an den Stellen (— 1,0) und (0,0) 


ER besitzt, ergibt sich der Wert des Integrals zu: 
Bild 4. 
=—-2ni|— Res. (0,0) — — Res.(-- 1, 0) 


Somit erhält man nach GI. (19), (20) und (21): 


für >0 


} 


In ganz entsprechender Weise läßt sich das Hakenintegral y (#, y; — 1) behandeln; man 
erhält schließlich das Ergebnis: 


2. yla,y; + _2niyetr'* (ve = 
v®- 


(22b). 


s-y 


VII. Die Geschwindigkeitskomponenten und das Potential. 


Die Werte der Gl. (22) für die Hakenintegrale können nun in die Gl. (19) für die 
Vertikalgeschwindigkeit v® eingesetzt werden. Man erkennt sofort, daß nur im Bereich 
x” —9y>0 die Geschwindigkeit » von Null verschieden ist; dieser Bereich stellt, nach der 
Definition von x und y (Gl. (17b)), ein unter dem Machschen Winkel 9 von der Flügel 
vorderkante ausgehendes Keilfeld dar. In diesem Bereich «—y>0 ergibt sich: 


10) W stsen: A treatise on the theory of Bessel funetions. Cambridge 1922, S. 359. 


v®-y 
| 
| 


| av 


lech. 
1942 


dt. 


oT. 


n) 


|). 
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x 


a 
o 


yeor* x le 


wobei in der oberen Grenze des Integrals über &, der obere bzw. untere Ausdruck zu nehmen 


ist, je nachdem ob # — 7-- y=0 ist, d. h. je nachdem der betrachtete Punkt (x, 9) im Be- 


2) z-y>0 2) 


&2-2, 


Bild 6. 


Bild 5. 


reich 1 oder 2 von Bild 5 liegt. Vertauscht man die Reihenfolge der Integrationen über & 
und x,, so ändern sich die Grenzen der Integrale in folgender Weise (s. Bild 6): 


1. y<0: \ \ 
2. 2-f—-y>0: \ \ 


Da uns zur Berechnung der Flügelkräfte nur der Bereich am Flügel selbst interessiert, 
führen wir die Berechnung nur in diesem Bereich 1 durch. Wenn man die Integrationen 
nach x, ausführt, die Konstanten a, , C,, C, nach Gl. (15) einsetzt, Gl. (16) berücksichtigt 
und auf das dimensionsbehaftete Koordinatensystem .r, y zurückgeht, so ergibt sich nach 
einiger Zwischenrechnung: 


z—t—-y<0 
ur} 
dr)x — yceotg* 
(23), 

u,sin d 
x | to? 9 

cos’ 


/(dg, 


y’cotg?’d 


wo und #=are sin =Machscher Winkel ist (24). 
o 


D, 

Um die Geschwindigkeitskomponente « zu bestimmen, müssen wir v nach . differen- 
zieren und nach y integrieren, siehe Gl. (5). Mit Hilfe der leicht zu verifizierenden 
Differentiationsrelation 


> 
x 
sy | 
| 
x 
£-y>0 
N 
7 | 
F | 
| 
), 
| 


.angew. Math. Mech. 


v. Borbely, Über die Luftkräfte bei Überschallgeschwindigkeit Nr.4 Aug. 19% 


a 73,2 
. 


S=by 


ergibt sich dann: 


u 


4, sind 


; 

ycotg 

11 


(m sin d 


- definiert ist, 


Für die Potentialfunktion 7, die dureh die Gleichungen 


erhält man: 
im Bereich 1 (Bild 5) 


x 


Ss=yeotg" 


cotg? 


u,sın d 


=yeotg* 


Damit sind die Geschwindigkeitskomponenten und das Potential im Bereich der oberen 
Halbebene (y >0) bekannt, jedoch noch nicht in der unteren, da die Hakenintegrale nur 
für positives y definiert sind. Man erhält jedoch die Werte in der unteren Halbebene 
leicht durch folgende Überlegung: Wenn längs der .x-Achse zur Zeit 7 bestimmte Vertikal- 
geschwindigkeiten herrschen, 
so haben sich diese nacheiner Y 
halben Schwingungsperiode, 
d. h. zur Zeit r+n/o auf 
den entgegengesetzt gleichen 
Wert geändert. Denkt man 
sich die ganzeobereHalbebene 
im Zustand r+7/o an der 
‚-Achse gespiegelt, so stim- 
men längs der .„-Achse die 

Vertikalgeschwindigkeiten 
mit denen der oberen Halb- e_ 
N 


ebene im Zustand r nach 
Größe und Richtung überein. 
Daraus folgt, daß der Strö- 
mungszustand der oberen 
Halbebene zur Zeit 
durch Spiegelung auf die 
untere Halbebene übertragen, 
die Strömung in der unteren v@;-yr) 
Halbebene zur Zeit r gibt 

(siehe auch Bild 7). an 
der Oberseite des Flügels 
also Störgeschwindigkeit und e 
Potential den Wert haben '*): Bild 7. 


11) An sich kann bei der Integration noch eine additive, willkürliche Funktion von x auftreten. In unserem 
Falle muß sie jedoch verschwinden; dies ergibt sich z. B. aus der Bedingung, daß beim Durehtritt der Strömung durch 
die Front der Machschen Welle, d. h. auf der Geraden x—yeotg®, die Geschwindigkeitskomponente in Richtung 
dieser Geraden ungeändert bleiben muß. 


12) Da die Schwingungsamplituden des Flügels klein sind, stimmt die Geschwindigkeit am Flügel in erster 
Näherung mit der an der x-Achse überein. 
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w=+u(e, 
v>+0 


so haben sie an der Unterseite den Wert: 
— y;t) 


v>—V 


IX. Kräfte und Momente am Flügel. 
Um die Flügelkräfte zu berechnen, muß man den Druck p am Flügel kennen. Wenn 
man in der Bernouillischen Gleichung: 


die durch Integration der Gl. (4) erhalten wird, die Zustandsgleichung p =«:o* einsetzt, so 
erhält man für p: 
oder, unter Beibehaltung nur der Glieder erster Ordnung: 
»—1o, 
Für die Auftriebskraft am Flügel ergibt sich dann: 
t t 
für das Moment der Luftkräfte um die Vorderkante (im mathematisch positiven Sinn ge- 


rechnet): 
t 


- 


x=0 


Setzt man hier « (2,0; 7) und we nach Gl. (25) und (26) ein und führt in die Formeln 
die momentane Auslenkung des Flügels y= Yei’" und den momentanen Anstellwinkel 
e=Eei’(—"o) ein (siehe Gl. (9)*), so erhält man 
1 
P 
= 


30, 


1 


0 
1 u 


\ 


u=0 


1 u 
\ e-iv:J, 
1 u 
\ u \ -\ \ (o sin 
u=0 u=0 


13) Die Glieder Eo, Yo der Gl. (9) verschwinden hier, wenn ınan zur Anströmrichtung symmetrische Schwin- 
gungen betrachtet. 


Mech. 
| 
26). 
ren 
nur 
ene 
kal- 
0 
. (da 
-:. 
ırch 
ung 
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oder schließlich nach einigen Umformungen 


Pr (31a). 
2; 
(Ha + Xu N, + > dr 
Ganz entsprechend wird 
M de de 
> 


Hierin sind die Größen r,, q4 Beiwerte, die von einer reduzierten Frequenz » = -- 


und dem Machschen Winkel 9 abhängen: 


fi (sind, o) = (osind-u)w-du, 4=0,1,2,3 


r, f, (sin (sin o) 


r,— f, (sin d,o) -f, (sin d, 4. — fı (Sind, ©) — f, (sin #, (32) 
r, (sin d,o) 2 f, (sin ©) 1=2f, (sind, —3f, (sind, ©) 
+ f, (sin 9, ©) + f, (sin 9, 


Um eine Vorstellung von den Unterschieden zwischen der inkompressiblen Strömung 
und der Überschallströmung zu geben, sind in Tafel 1 die Luftkräfte beider Strömungen 
gegenübergestellt. In der gewählten Darstellung haben bei inkompressibler Strömung die 
Kräfte feste Angriffspunkte, während sie bei Überschallströmung im Lauf einer Schwingung 
ihre Lage ändern (die Hebelarme sind komplex); ferner ist im ersten Fall die Bewegung 
verschiedener Flügelpunkte (vorderer und hinterer Neutralpunkt, Mittelpunkt) für die 
Kräfte maßgebend, während sich im zweiten Fall alles auf die Bewegung der Vorderkante 
beziehen läßt. 


X. Berechnung der Beiwerte r, y. 

Für die Integrale f; der Gl. (32), die zur Ermittlung der Werte r, 4 numerisch berechnet 
werden müssen, läßt.sich eine Rekursionsformel aufstellen, die f,, f;, f, auf f, zurückführt. 
Führen wir den Ausdruck ein 


F; (sin 9, 


(u-sind)du. . 2 
0 


so erhält man durch partielle Integration leicht: 


F, (sn J, (u sin d) iA F,_Msind,o)+isind (usind)du 


und 


\ein (usind) du =ie-i"wJ, (usin ®) isin F, (sin o) 


(usind)du. 


Führen wir die letzte Gleichung in die vorletzte ein und benutzen die vorletzte nach Ver- 
tauschung von 4 mit A — 1, so erhalten wir die Rekursionsformel: 


(1 sin? 9) (sin ei” J, sin d) — sin @4.J, sin 
(0 sind) +(1— 24) F,_, (sin 9, ©) ei® (34). 
+(1-— F,_, (sin 9, ei® , i>1 


Für =1 ergibt sich speziell 
(1 sin’) F, (sin d, —=i oJ, sind) - »sindJ, (osin®) iF,(sin d,o)eiv , (34a). 


Damit lassen sich aus dem Integral F, bzw. f, alle F} bzw. f} berechnen. 


ot 
XUM 
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Gegenüberstellung der Luftkräfte bei inkompressibler Strömung und Überschallströmung. 


Inkompressible Strömung 
(nach Kassner-Fingado, |. «e.) 


Überschallströmung 


1. Eine dem wirksamen Anstellwinkel 
‚les hinteren Neutralpunktes propor- 
tionale Auftriebskraft 


P = Auf: 
triebs- 
vektor 

(siehel. ce.) 


= ‚P 


Angriffspunkt der Kraft im vorderen 
Neutralpunkt: 


1. Eine dem wirksamen Anstell- 
winkel der Vorderkante pro- 
portionale Auftriebskraft 


| 1 dy. 
Angriffspunkt der Kraft: 


r, 


2. Eine der Drehgeschwindigkeit des Flügels pro- 
portionale Auftriebskraft 
K de 
o 
dr - 


Angriffspunkt der Kraft im hinteren Neutral- 
punkt: 


,=0,75t. 


2. Eine der Drehgeschwindigkeit 


des Flügels proportionale Auf- 
triebskraft 
de 


K,= 4 "dr 


Angriffspunkt der Kraft: 


Ya 


1 
2 


. Eine der Beschleunigung des Flügelmittelpunktes 
proportionale Auftriebskraft 


0, 8? 


Angriffspunkt der Kraft im Flügelmittelpunkt: 


.Eine der Beschleunigung der 
Vorderkante proportionale 
Auftriebskraft 


„4? 


RK, = 4 de: -P; 


Angriffspunkt der Kraft: 


h,= t. 


Ta 


4.Ein der Drehbeschleunigung des Flügels pro- 
portionales Moment 


d?e 
- 


. Eine der Drehbeschleunigung 
des Flügels  proportionale 
Auftriebskraft 


Angriffspunkt der Kraft: 


19 
Br; 


Die Berechnung von f, wird mit Hilfe von Reihenentwicklung vorgenommen. 


wir ın - 


Führen 


f, (sin = = \ e-'uJ, (usind)du 


0 


die Produktenentwicklung der Besselfunktionen 


fein _ Dil 


ein, so erhalten wir 


1 
f, (sin d, ©) = = > dl 


14) Siehe z. B. Watson, |. ce. S. 142, Gl. (16). 


— sin? 


Al 


etw-u) J,(u)udu. 


0 
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34). 
| 
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Das in dieser Reihe auftretende Integral läßt sich ausrechnen. Durch partielle Integration 


erhält man 


J, (u); 


Aryl 
0 


lau 


\ J,(u)du 


art 


ir 


oder 


(244 1) \ \ui+t (wdu. 
0 


Integriert man das zweite Glied der rechten Seite noch einmal partiell: 


so erhält man schließlich 


[07 


0 


Damit ergibt sich 


i=0 


Diese Formel läßt sich zur numerischen Berechnung verwenden. 


XI. Anhang: Die stationäre Strömung. 


Der Vollständigkeit halber soll noch gezeigt werden, daß sich auch die stationäre 
Strömung um ein Flügelprofil, die an sich schon bekannt ist'*), mit Hilfe von Hakenintegralen 
behandeln läßt. 


Wir bauen das Strömungsfeld wieder aus Elementarströmungen auf. Diese ergeben sich 
aus Gl. (11), indem wir dort 6>0 gehen lassen, d. h. die Strömung stätionär werden lassen: 


U. 


u 


Wie sich leicht zeigen läßt°), hat v, je nach der Lage des Punktes x, y die Werte 


0 für — 


® 
| 1 für ee, — yootgd >0 


d. h. », hat in der «-y-Ebene nur in einem vom Punkt x, ausgehenden Keil (mit dem halben 
Öffnungswinkel d) einen von Null verschiedenen, überall konstanten Wert. Das von der 
Elementarlösung dargestellte Strömungsfeld entspricht also der Prandtlschen Strömung um 
eine Ecke (bei sehr kleinem Neigungswinkel der Ecke). 


18). Siehe z. B. Ackeret: Luftkräfte auf Flügel, die mit größerer als Schallgeschwindigkeit bewegt werden. 
Z. Flugtechn. Motorl. Bd. 16 (1925), S.72. Prandt]: Allgemeine Überlegungen über die Strömung zusammendrückbarer 
Flüssigkeiten. Volta-Tagung, Rom 1936. 
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Wir überlagern nun derartige Elementarlösungen, so daß sie die Randbedingungen der 
Strömung um einen Flügel erfüllen. Diese lauten (s. Abschn. V und Gl. (10b)): 


längs der «-Achse 1) für 
2) für U, E, + Ö,, u 
3) für t<. v=0. 


Folgender Ansatz führt zum Ziel: 


1 ei HIX— X — yeotg‘') 
\ du \ 


ei yeotg®) 
(x,y) = uld., 


a D, ein (x— yeotg®) eiu(x —t— yeotg®) 

0 


Wie man leicht sieht, sind die Randbedingungen 1) und 3) ohne weiteres erfüllt; damit 
auch Bedingung 2) erfüllt ist, müssen die noch freien Größen f («,) und € folgender Gleichung 
genügen: 


x 


\ f(=,) de, + 


d? > d Ö,, u 
U, (@); ( =— dr AR 


Der Index o,u bei ö soll andeuten, daß im oberen Teil der Halbebene ,, im unteren Teil ö, 
in Frage kommt. Es ergibt sich dann: 


im Bereich 0< — |y| cotgd<1t 


E, 


Für die andere Geschwindigkeitskomponente « erhält man nach Gl. (5): 


dö 

wobei das negative Zeichen für die obere, das positive für die untere Halbebene gilt. 
Die Berechnung der Flügelkräfte geht dann in bekannter Weise vor sich. Es sei nur 

kurz das Ergebnis angegeben, das sich genau mit den Formeln von Ackeret deckt: 


Auftrieb 


Moment um die Vorderkante 


(im math. posit. Sinn) — 


—=20,03 tg9d- 


Widerstand WER-P: 


Abgeschlossen im Oktober 1939. 
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Über eine gesteuerte eindimensionale Bewegung. 
Von Herbert Bilharz in Braunschweig. 
Bericht aus dem Institut für Aerodynamik der Luftfahrtforschungsanstalt Hermann Göring. Braunschweig'*). 


Die periodischen Lösungen der nichtlinearen gewöhnlichen Differential- 
gleichung (1) werden explizit angegeben und anschließend die Stabilitäts- 
verhältnisse der geradlinigen Bewegung eines Massenpunktes bei unstetiger 
rerallgemeinerter Kraft (3) untersucht. Solche Bewegungen treten auf bei 
Benutzung bestimmter fehlerfreier Automatiken zur Stabilisierung dyna- 
mischer Systeme eines Freiheitsgrades. Wie Verfasser an anderer Stelle 
für u—=0 in (1) gezeigt hat. können fehlerbehaftete Automatiken die Ergeb- 
nisse u. U. verbessern. 


I. Einleitung. 
Es sei .(f|o) eine reelle, zweimal abteilungsweise stetig differenzierbare Funktion der 
unabhängigen Variablen £>0 und des reellen Parameters o mit -w<o<x, 


ir’? >0,2>0 und « >0 seien Konstante. Unter diesen Voraussetzungen werden im 


folgenden die Integrale der nichtlinearen gewöhnlichen Differentialgleichung zweiter Ordnung 


untersucht; die Funktion sgn + 0x) ist definiert durch 


1 0 


(sgn 0 wird also nicht erklärt. Für «=0 habe ich anläßlich einer Stabilisierungsaufgabe 
der Flugmechanik ') versucht, allgemeine Aussagen über die Lösungen von (1) zu erhalten. 

Deuten wir f als Zeitvariable, dann ist (1) die Differentialgleichung für die geradlinige 
Bewegung des Massenpunktes Eins in einem verallgemeinerten Kraftfeld mit der unstetigen, 
von der Ortskoordinate x, von der Geschwindigkeit # und dem Parameter o abhängigen Kraft 


Wird darin y=. als neue Variable eingeführt (und „7 für 7 geschrieben), so ist (1) äqui- 
valent dem System von zwei gewöhnlichen Differentialgleichungen erster Ordnung 


Nach den geometrischen Methoden Poincares?) lassen sich alle möglichen Bewegungen 
solcher dynamischen Systeme mit einem Freiheitsgrad qualitativ untersuchen; den periodischen 
Lösungen (mit Einschluß der Gleichgewichtslagen) kommt darnach wesentliche Bedeutung 
für Stabilitätsaussagen zu. Da Gleichgewichtslagen für (1) nicht definiert sind, suchen wir 
periodische Bewegungen aufzufinden und werden sie im nächsten Abschnitt explizit angeben. 
Damit wird insbesondere die Existenz von Lösungen der Gl. (1) für alle £>0 bewiesen sein. 
Für niechtlineare Differentialgleichungen ist diese explizite Angabe der periodischen Lösungen 
im allgemeinen nicht in geschlossener Form möglich®). Hier, wie bei einer großen Anzahl 
ähnlicher Aufgaben, ist sie besonders einfach; im wesentlichen ist nur ein lineares Gleichungs- 
system aufzulösen. 


*), Im Auszug vorgetragen auf der GAMM-Tagung 1942 in München. 


!) Luftf.-Forschg. Bd. 18 (1941), S. 317 bis 326; im folgenden zitiert mit [L]. Dort sind auch einige 
technische Hinweise gegeben. Vgl. ferner K.Bögel: Das Verhalten gedämpfter und aufschaukeln- 
der freier Schwinger unter der gleichzeitigen Einwirkung einer konstanten Reibungskraft. Ing.- 
Arch. Bd. 12 (1941), S. 247 bis 254. Darin wird in (1) unsere Signumfunktion durch die speziellere 
Funktion sgn # ersetzt. Die Ergebnisse von Bögel lassen sich aus unseren Resultaten ableiten, 
wenn wir 0 —=x setzen; die dortige kritische Schwingungsweite entspricht unserer periodischen 
Lösung. 

H. Poincare&: Memoire sur les courbes döfinies par une &quation differentielle. Journal de Mathe- 
matiques, ser. 3, vol. 7 (1881), vol. 8 (1882), ser. 4, vol. 1 (1885), vol. 2 (1886): insbesondere die 
beiden erstgenannten Teile. 

Vgl. z. B. G. Hamel: Erzwungene Schwingungen endlicher Amplituden. Math. Ann. Bd. 86 (1922), 
S. 1 bis 13. 


IE”). 


tial- 
tät 
iger 
f bei 
telle 
"geb- 


(2), 


zabe 
Iten. 
nige 
gen, 
‚raft 


iqui- 


(4). 


Igen 
tung 

wir 
ben. 
sein. 
zahl 
ngs- 


inige 
keln- 
Ing.- 
llere 
iten, 
chen 


ıthe- 


Au Nr.4 Aug. 1942 Bilharz, Über eine gesteuerte eindimensionale Bewegung 207 


Mit Benutzung einer Stabilitätsdefinition von Trefftz wird im letzten Abschnitt ge- 
zeigt, daß die periodischen Lösungen von (1) (für jedes o mit -o<o<0 gibt es genau eine 
und für o>0 keine) stabil sind und alle anderen stabilen Bewegungen für >» gegen diese 
periodische Bewegung konvergieren. Instabile Lösungen treten dadurch auf, daß von einer 
endlichen Zeit #, an (1) nicht eindeutig fortsetzbar wird‘). Nach Befragen des Experiments 
und Berücksichtigung seiner Fehler ((2) ist praktisch nur näherungsweise zu verwirklichen) 
läßt sich eventuell eine eindeutige Fortsetzung für #>t, angeben. Derartige Versuche für 
«a >0 sind mir bis Jetzt nicht bekannt geworden; ich muß mich darauf beschränken, auf diese 
Instabilität hinzuweisen. 

Für die Integrale der Systeme 


bzw. 


mit 


und reellem /" gelten die bekannten Existenz- und Eindeutigkeitssätze°), von denen später, 
insbesondere beim Stabilitätsnachweis, öfters Gebrauch gemacht wird, ohne daß diese Tat- 
sache jedesmal ausführlich Erwähnung findet. 

Aus der Schwingungslehre kommen die hinreichend bekannten Tatsachen zur Anwendung: 
Die Integrale & (f), y(f) dieser Systeme geben in der (x, y)-Ebene (Phasenraum) orientierte 
doppelpunktfreie, diese Ebene einfach und lückenlos überdeekende Kurven, welche sich für 


I 
#»’< u abzählbar unendlich oft um die Punkte = bzw. % (Strudelpunkte) herumwinden und 


für in diese konvergieren. Die Richtung der Phasenkurven ist die »-Achse 
wird also stets vertikal überschritten. Sei Pf) ein laufender Punkt längs einer solchen 
Spirale, so ist der Abstand zwischen P und dem Strudelpunkt eine monoton gegen Null 
konvergierende Funktion von # für >». Für #° > u münden die Phasenkurven für t>x 


ebenfalls in die (Knoten-) Punkte „ ein, überqueren die «-Achse aber höchstens einmal. 


Nachstehende Definition ist nicht unentbehrlich, aber oft zweckmäßig: Die Gesamtheit 
aller Integrale x(t|o) von (1) für festes o heißt ein gesteuertes System I,, das einzelne 
Integral die Lösung in einem gesteuerten System. 


II. Periodische Lösungen. 

Nach (2) ist y-sgn(e-+ox&) eine abteilungsweise konstante Funktion von ? (und natür- 
lich auch von o, worauf wir später zurückkommen werden); die Länge der einzelnen Ab- 
teilungen wird durch & (t|o) bestimmt. Periodische Lösungen von (1) werden sich dort ein- 
stellen, wo diese Abteilungslängen konstant sind. 

Wir bestimmen nunmehr eine Klasse periodischer Lösungen und wählen hernach die- 
jenigen von (1) aus. Es sei a >0, s durchlaufe die Menge der ganzen Zahlen ..., —2, - 10, 
+1, +2, ....; ferner sei /'(f) eine reelle, abteilungsweise konstante und mod. 2a periodische 
Funktion von mit |7'|—=y>0. Speziell werden gesetzt) 


2sa<t<(2s+l)a, 


Die Zeitmarke t=0 ist danach ohne Beschränkung der Allgemeinheit eindeutig in einer 
Sprungstelle von — /' nach /' (kurz: Aufschaltstelle von /') festgelegt. Diese Normierung 
bleibt vorläufig bestehen. Die t-Intervalle rechter Seite, in (8) seien mit J,, bzw. J,, be- 
zeichnet; J,,—+ J,, ist das Grundintervall. 

An Stelle von (1) betrachten wir die gewöhnliche lineare, inhomogene Differential- 
gleichung zweiter Ordnung mit konstanten Koeffizienten 


*) Über das Auftreten solcher Stellen vgl. die Bemerkungen bei Painlev&: Lecons sur la theorie 
des &quations differentielles. Paris 1897, S. 18 der Introduction. 

5) Vgl. z. B. E. Goursat: Cours d’Analyse. Bd. 2, 5. Aufl. Paris 1933, S. 364 bis 384. 

) Mit!=t Ist (8) durch die für alle konvergente Fouriersche Reihe +1 
darstellbar. Hiervon sehen wir bewußt ab. ver 
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und suchen hierfür (periodische) Lösungen. Alle Schaltstellen von /' sind auch solche von 
y-sgn(c+0x%), aber nicht umgekehrt; (1) und (9) sind also im allgemeinen nicht äquivalent. 

Sind /,=+4, die beiden reellen oder konjugiert komplexen Wurzeln der charakteristischen 
Gleichung 


so bilden bekanntlich eı* und e‘e! ein Fundamentalsystem der homogenen Gleichung 
Für = u, also A), =A,—= — x, sind und t-e”*! zwei linear unabhängige Integrale von (11). 


Partikularintegrale von (9) sind 2=— im J, und J,s. Mit den zunächst unbe- 


stimmten Konstanten 4A,,, B,s, A,s und B,, haben darnach die allgemeinsten Lösungen von (9) 
die Form 


und 


Zur Bestimmung dieser Konstanten werde gefordert: 


1. Stetigkeit der Integrale (12) und ihrer ersten nach f genommenen Ableitungen y(f) 
in den Intervallen J,;+J,;, also 


y(2s+Da-0)=y(2s+Da+0) (18); 
2. Periodizität mod. 2a in ? von x und 9, also 
—0) 


In jedem der Fälle #°= « ergibt sich hieraus je ein lineares Gleichungssystem für die vier 
Unbekannten A,,, ..., B,s mit den Koeffizientendeterminanten 


(14). 


— A,)? e-2artas +1) +3) — — +1) für . (15) 
und 


Beide Determinanten sind reell und von Null verschieden, D, wegen a >0, A, =4, und D, 
wegen des reellen ax=F0. Die Konstanten werden 


Bs = 27 (A, — A,) (ea +42) — — (16,) 
s D; u 1 
Aserch, 
für u und entsprechend 
DER. e?uns 
Au= 
1) 
(16,) 
für u. 
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Mit den Abkürzungen 


für 
T=+yr—u 
folgen schließlich aus (16) nach einiger /wischenrechnung, für Anfangslage «,—=x(0) und An- 


fangsgeschwindigkeit y,— y(0) als Funktionen von a die in nachstehender Tabelle zusammen- 
gestellten Ausdrücke 


(17) 


I zsinao—oZinax : — Sinax »Sinar— ı&inax 
uo 1+Coax ur 


o cosac-+ Gojax 1+Cojax Cojar+bojax 


(>00) 


Die in Klammern beigefügten Zahlen geben den entsprechenden Limes an für a>x. Wird 7' 
durch — /' ersetzt, so ändern x, und ,, nur ihre Vorzeichen. 

‚ Im Intervall J,, erhalten wir mit den Anfangswerten (18) für die Integrale « (f) von (9) 
und deren Ableitungen (f) 


yh=e-t In cos ot sin 


für #’< u; 


r 


r 
le, + +5. | 


r 
| 


+ 
Yo 


für #’— u; 


für > u. 

Diese Lösungen lassen sich unter Zuhilfenahme von (16) für £ > «a stetig fortsetzen und 
genügen für ganzzahliges n den Funktionalgleichungen 


| 
) 
haben also die Periode 2a und den Mittelwert Null. 
Zu jedem a >0 gibt es ein reelles o mit 
Wegen (20) bestehen für ganzzahliges n>1 dann auch die Gleichungen 
z(na)+oy(in a)—=V 
Aus (21) und (18) folgen die von /’ unabhängigen Ausdrücke 
1 Sina 
sinao, 


a» Sinax 


Sinaz 
omar 


| 
un 
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Bild 1. 


Bild 1 zeigt für 7’>0 in jedem der Fälle #°Z u den Verlauf der Funktionen «,, y, aus 
(18) und o aus (23) in ihrer Abhängigkeit von a >0, der einzigen noch Unbestimmten. Dieser 
Darstellung ist zu entnehmen, daß für #° > « die Ergebnisse sich nur quantitativ von den- 
jenigen für z#°— u unterscheiden. 

Die Funktion 

(24) 
hat nach (22) die Eigenschaft, an ihren (reellen) Nullstellen t=n.«a den Übergang von T 
nach — /' und umgekehrt zu steuern. An den Stellen »a werden also die beiden Funk- 
tionen f; aus (7) ineinander übergeführt. Wenn {(t) nur die Nullstellen na hat (also die 
Intervalle na<t<(n+ 1)a nullstellenfrei sind), stimmen für /'—=y die periodischen Lösungen 
von (1) mit denen von (9) überein. Wir werden später sehen, daß zwischen den Wurzeln 
na und (n-+1)a der Gleichung {({l)=0 noch weitere Wurzeln liegen können. 

Vermöge der Lösungen (12) lassen sich an- 
schließend die Fragen nach der Existenz ge- 
steuerter Systeme 2o, in welchen es periodische 2 
Lösungen gibt, nach den Anfangswerten dieser 
Lösungen und den zugehörigen Perioden mühe- 
los beantworten. Die hierzu notwendigen Auf- 
lösungen nach a der transzendenten GIn. (18) 
führen wir graphisch aus. In Bild 2 sind in der 
(2, , Ebene (kurz: Nullebene) die, Kurven (%,) 
und a(y,) aus Bild 1 wiedergegeben. 

Für #’< u sind diese Funktionen unend- 
lich vieldeutig. «x, (y,) ist eine den Nullpunkt mit 
vertikaler Tangente verlassende Spirale, die sich 


abzählbar unendlich oft um den Punkt (= .o) 


herumwindet und schließlich in jenen Punkt Bild 2. 
konvergiert. Die x,-Achse wird stets vertikal 
überschritten. Die Halbperiode a (y,), die auch auf der #,-Achse aufgetragen ist, ergibt eine 
sich bis ins Unendliche erstreckende Kurve in Gestalt einer gedämpften Welle, welche die 
&,-Achse ebenfalls abzählbar unendlich oft schneidet. Unter Benutzung nachstehender Vor- 
schrift können die Punkte der Kurven x, und a umkehrbar eindeutig einander zugeordnet 
werden: Vom Nullpunkt ausgehend und auf der &,- bzw. a-Kurve fortschreitend, numeriere 
man, jedesmal mit 1 beginnend, die Schnittpunkte dieser Kurven mit der «,-Achse. Ferner 
ordne man die Bogenstücke dieser zwei Kurven zwischen je zwei gleichnumerierten Paaren 
konsekutiver Schnittpunkte einander zu, wobei die Bogenlängen vom Schnittpunkt kleinerer 
ÖOrdnungszahl zu demjenigen größerer Ordnungszahl positiv gezählt werden. 

Für #’> u sind «, und a zweiwertige Funktionen von 9,: x, und a lassen sich aber 
auch hier nach obiger Vorschrift eineindeutig einander zuordnen. 

Gl. (21) repräsentiert in der Nullebene ein durch den Koordinatenursprung gehendes 


Geradenbüschel; die einzelne Gerade @ dieses Büschels (festes o) hat die Steigung — S 
Neigungswinkel ö bestimmt sich eindeutig aus der Gleichung 
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Die vom Nullpunkt verschiedenen Schnittpunkte von @ mit der x,-Kurve geben Anfangs- 
lagen und nach Bild 2 Anfangsgeschwindigkeiten und Intervallängen der Lösungen (12). 
Für #°<u« erhalten wir bei Beschränkung auf endliches o solche Schnittpunkte für 


« <P ist dabei die Neigung der Tangente aus dem Ursprung an den dem Intervall 
a<ao<2n zugeordneten Bogen der Kurve x,. Je nach der Wahl von o im Intervall (26) 
gibt es eine oder zwei oder drei usw., jedenfalls abzählbar viele Lösungen und keine für 
0<o< 0o*. 


Ist #°> u, dann gibt es je eine Lösung für 
> 


und keine für <o<x. 

Aus diesen Lösungen von ©) sind nunmehr die periodischen Lösungen von (1) aus- 
zusondern; im folgenden bleibt es freilich dahingestellt, ob den Lösungen mit sehr großer 
Periode irgendeine physikalisch-technische Bedeutung zukommt. 

Bislang war t=0 Aufschaltstelle von 7. Die Zeitvariable werde (durch evt. Addition 
einer Konstanten) so normiert, daß t=0 Aufschaltstelle von y(>0) ist. Diese Normierung 
wird von jetzt ab stets beibehalten. 

Es sei # >0 beliebig klein. Eine Lösung von (1) mit den Anfangswerten x,, y, existiert 
nur dann für O<t<e (d. h. ist über {=0 hinaus fortsetzbar), falls die Steuerungsfunktion 
(24) mit Z(0)=0 für #>0 zunächst abnimmt. Es muß also die Ungleichung 


t=+V0 
bestehen, daraus folgt 


> 


In der Nullebene gibt Z,=0 einen Kegel- 
schnitt (Ellipse, Parabel bzw. Hyperbel für 
»Zu), der die y, Achse im Ursprung tan- 


giert und durch den Punkt (2, 0) hindurch- 


geht. In Bild 3 sind diese Kurven dargestellt; 
die dort schraffierten Gebiete geben den Werte- 
vorrat von (28) an’). Durch diese Anfangs- 
werte gehen Lösungskurven von (1), die vor- 
erst nur in einem kleinen Zeitintervall defi- 
niert sind, während es uns darauf ankommt, 
Lösungen zu erhalten, deren Existenz für alle 
t>0 gesichert ist. 

Übertragen wir die «,-Kurven von (9) 
aus Bild 2 in Bild 3, so scheiden die Anfangs- 
werte x,, 9% von (9) mit "=y>0 für 
*<o<mw in (26) aus, weil sie die Un- 
gleichung (28) nicht erfüllen‘). Die verbleibenden Punkte der «,-Kurven liegen sämtlich im 
Quadranten x,<0, y,<_ und genügen (28); Anfangswerte der periodischen Lösungen von (1) 
sind sie jedoch nur dann, wenn Z(f)<0 für O<t<a ist. 


Bild». 


‘) Für #? > u hat die Hyperbel £,—=0 einen ganz im zweiten Quadranten gelegenen Zweig, der vom 
Nullpunkt endlichen Abstand hat; im Innern dieses Zweiges ist (28) erfüllt. 

Daß (1) für 0E&>0 und #°< u keine periodische Lösung mit a >0 hat, läßt sich nach einer münd- 
lichen Mitteilung von H. Seifert auch folgendermaßen einsehen: Der Nullpunkt U der Phasen- 
ebene halbiert bei einer periodischen Lösung nach (20) die Strecke der Steuerungsgeraden zwischen 
!=0 und t=a, wobei der durch Steuerungsgerade und Phasenkurve in diesem Zeitintervall be- 
grenzte Bereich konvex ist. Für die Aufschaltstelle von 7 >0 liegt U links vom Strudelpunkt S. 
Wegen der monotonen Abnahme des Abstandes eines laufenden Punktes P (£) der Phasenkurve vom 
Strudelpunkt ist SP(a)< SPD), also der Winkel SU P (a) kleiner = . Da US die positive z-Richtung 


angibt (die positive Y-Richtung geht daraus durch Drehung um x/2 um den Ursprung U hervor), 


liegt die Steuerungsgerade mit der Steigung — = notwendig im ersten und vierten Quadranten, 
d.h. es ist o<0. 


x 
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Wegen £(0)=L(a)=0 genügt es zu zeigen, daß £(f) im Intervall O<t<a nur eine 
Nullstelle besitzt und die Bedingung (28) erfüllt. Durch direkte Rechnung aus (18), (19) und 
(23) folgt im Fall u, daß monoton steigend verläuft für jedes a >0. Ist <u, so 
gilt zunächst 


Die zwischen x und 27 gelegene obere Schranke 4 für ao ergibt sich aus dem Schnittpunkt 
der Ellipse 


mit der x,-Spirale in der Nullebene. (31) folgt aus (29) durch Vorzeichenumkehr von «, 
und %,. Die Phasenkurve (s. u.) der (periodischen) Lösung von (9), deren Anfangswerte durch 
diesen Schnittpunkt definiert sind, tangiert nämlich R 
—1 
die zugehörige Steuerungsgerade der Steigung — 
mit o>0o*. Für ao>4A haben die Phasenkurven 
Schleifen, d. h. £(t) hat im Intervall O<t<a 
mindestens eine Nullstelle.e. Um außer (30) noch 
die Ungleichung (28) zu erfüllen, muß schließlich 
die Ungleichung 


bestehen. In Verbindung mit Bild 1 ist also nur 
der erste äußere Teilbogen der x,-Kurve geometri- 
scher Ort der Anfangswerte periodischer Lösungen in den gesteuerten Systemen Io. Wegen 
y>Pb ist -—o<o<0. Bild 4 zeigt die Phasenkurve einer ausgeschiedenen Lösung von (9) 
mit 27<ao< 37; sie ist keine periodische Lösung von (1). 


Bild 4. 


Wir fassen die bisherigen Ergebnisse zusammen: In den gesteuerten Systemen Io mit 
— o<o<O hat (1) genau eine periodische Lösung mit explizit angebbaren Anfangswerten: 
für #’< u ist mit o aus (17) die Halbperiode a durch (32) beschränkt. 


Die weitere Bedeutung dieser speziellen für alle ?>0 definierten Lösungen von (1) er- 
hellt aus dem nächsten Abschnitt. 


II. Stabilitätsnach weis. 


Neben (1) betrachten wir dds äquivalente System (4) im Phasenraum°). Dieser ist 
hierfür eine Ebene '’), wenn als Koordinaten in einem rechtwinkligen kartesischen Achsen- 
system Lage x und Geschwindigkeit y aufgetragen werden. x (f) und y(t) ist die Parameter- 
darstellung der mit f orientierten Phasenkurve, ihre Neigung bestimmt sich durch 


dy_f 


(33). 


Die durch den Koordinatenursprung gehende Steuerungsgerade Z(t)=V der Steigung — - 


teilt den Phasenraum in zwei Halbebenen 
. 


In diesen hat nach (2) die verallgemeinerte Kraft (3) die Werte f+ aus (7), worin nunmehr y 
an Stelle von I’ zu schreiben ist. 


Das Integral x (t|o) von (1) ist eine abteilungsweise zweimal stetig differenzierbare 
Funktion, die aus den Lösungen der Systeme (5) und (6) auf der Steuerungsgeraden zusammen- 
gesetzt ist; die Ableitungen in Punkten dieser Geraden existieren einseitig. Der Nullpunkt 
=y=V0 ist auszuschließen, weil für ihn die Signumbedingung nicht definiert wurde. Für 
eine stabile Bewegung x (f|o) werden wir vernünftigerweise fordern, daß ihre Phasenkurve 


)YV gl. G. D. Birkhoff: Einige Probleme der Dynamik. Jhber. Deutsch. Math.-Vereinig. Bd. 38 (1929), 
S. 1 bis 16. 
'°) Diese Ebene ist doppelt überdeckt und ähnlich einer zweiblättrigen Riemannschen Fläche längs 
der Geraden £(f)=0 durch Null und Unendlich verheftet zu denken. 


| 
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für alle # definiert sei und ganz im Endlichen verlaufe, also z. B. ganz im Innern des Kreis- 
ringgebietes 

liege bei beliebig großem, jedoch endlichem R >0. 

Die Phasenkurven der im Abschnitt II angegebenen periodischen Lösungen sind end- 
liche, abteilungsweise stetig differenzierbare, geschlossene doppelpunktfreie Kurven; sie ge- 
nügen sämtlich der notwendigen Bedingung (35). 

Nach Voraussetzung ist die unabhängige Variable £ monoton wachsend, die Bewegung 
also progressiv; über die retrograden Bewegungen, bei denen —? monoton wächst, seien 
einige Bemerkungen kurz eingefügt: 

Der Übergang von f nach —t in (1) führt auf die Gleichung 


Andererseits folgt (36) aus (1) durch die Substitution 


. 


Dieses Resultat bestätigen auch die Formeln (18) und (23), in (18) ist aber zu beachten, daß 
die Zeitmarke t=0 der retrograden Bewegung in einer Abschaltstelle von y liegt. Unter 
Zuhilfenahme dieser Bewegungen lassen sich die Lösungen von (1) diskutieren, wenn die 
ursprüngliche Voraussetzung x >0 (die auch ferner beibehalten werden soll), durch <<0 er- 
setzt wird. Da die unabhängige Variable in (1) nicht explizit vorkommt, ist im Phasenraum 
die Orientierung der retrograden Bewegungen derjenigen der progressiven Bewegungen: ent- 
gegengesetzt. 

Abweichend von [I] legen wir für das Folgende eine von Trefftz'') gegebene Sta- 
bilitätsdefinition zugrunde: Eine Lösung der Differentialgleichung (1) heißt stabil, wenn sie 
für alle t>0 erklärt und mit Einschluß ihrem ersten Ableitung nach t beschränkt ist, also 
(35) erfüllt, und jede hinreichend kleine, aber endliche Störung dieser Bewegung abklingt. 
Diese Definition ist den Bewegungen in nichtkonservativen dynamischen Systemen angepaßt; 
ein derartiges System liegt wegen #==0 in (1) vor. 

Eine hiernach stabile Lösung tritt niemals isoliert auf'?); falls die periodischen Lösungen 
von (1) stabil sind, müssen darnach in ihrer Nachbarschaft weitere stabile Lösungen vor- 
kommen. Dies ist in der Tat der Fall. Es gilt nämlich der Satz: Die in-Il. angegebenen 
periodischen Lösungen von (1) mit den Anfangswerten &,, 9, Sind stabil und die Nachbar- 
lösungen X (f) mit den Anfangswerten X,, Y,=F 9, die durch 

und 
mit beliebig großem R>0 und einer geeignet gewählten positiven Konstanten c=c (x, u, y) 
gegeben sein sollen, sind ebenfalls stabil und konvergieren für t>x gegen die zum ge- 
steuerten System I. gehörige periodische Lösung. 

Gl. (38) ist eine Folge unserer Zeitnormierung, (39) in Verbindung mit (38) enthält im 
wesentlichen das Trefftzsche Stabilitätsmaß. Die Anfangswerte X,, Y, liegen nach (38) 
auf der Steuerungsgeraden, aber nicht notwendig in infinitesimaler Nachbarschaft von &,, %,: 
Wegen (27), (26), (32) und (39) ist -o<o<V0. 

Die Phasenkurve der periodischen Lö- uf) 
sung des gesteuerten Systems 2, durch die 
Anfangswerte x,, , sei die Kurve ABCDA 
in Bild 5. Für t>0 tritt die Phasenkurve 
der Nachbarbewegung X (t) durch die Anfangs- 
werte P,(X,, Y,) in die Halbebene £ (f) <0 und 
genügt dem System (6). Wegen Y, < 0 nimmt 
X(f) zunächst ab, erreicht für Y(t)=0 ein end- 
liches Minimum >— Rin P,(X<0),Y=0); 
anschließend verläuft die Phasenkurve in der 
Halbebene Y>0, X(f) nimmt zu. Zur Zeit 
t=t,(>0) wird die Steuerungsgerade in 
P,(X >0, Y(t,) >0) erreicht, um schlie£- 


a, E. Trefftz: Zu den Grundlagen der Schwingungstheorie. Math. Ann. Bd. 95 (1925), S. 307 bis 312. 
Dort auch die analytische Formulierung obiger Definition. 
'?) Für den Beweis vgl. Trefftz in loe. eit.''), S. 309. 
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lich für *>x in den Punkt Y=0, = zu münden (Gleichgewichtslage von (6)). Dieser 
letzte Teil der Lösung (6) interessiert hierbei nicht. Wesentlich ist nur, daß P, und P, auf 
der Steuerungsgeraden durch den Nullpunkt X=0, Y=0 getrennt werden. Die Koordinaten 
von P, sind kleiner oder größer als die Absolutwerte der Koordinaten von P,, je nachdem 
ob der Bogen P,P, P, außerhalb oder innerhalb des Bogens ABC der periodischen Lösung 
verläuft (d.h. Y,=9,). Der Beweis hierfür folgt rein analytisch aus (19) und (23). Analog 
[Z, pp. 321-322] läßt sich nämlich die Lösung X (t) einschließen zwischen der periodischen 
Lösung in 2, und derjenigen eines gesteuerten Systems I,, mit der Anfangsgeschwindigkeit Y,; 
man beachte dabei die eindeutige Bestimmtheit von (6) in der ganzen Phasenebene, Wegen 
Y„=y, werden 0,=o und nach (33) a,=«a. Daraus folgen außer der Existenz einer Null- 


o— 


stelle von EX die Schranken 


Die eindeutige Fortsetzung der Lösung (1) für >#, wird durch das System (5) in die 
Halbebene {(f) >UV gegeben. Wir erhalten den Bogen P,P, P,, der zur Zeit t=t#,(>t,) die 


Steuerungsgerade trifft (evtl. nur einmalig für x’ = u) und für >» gegen | = 2 0) konvergiert. 


Für t >, interessiert diese Bewegung nicht. Wird vergleichsweise die periodische Lösung 
eines gesteuerten Systems Io, mit der Anfangsgeschwindigkeit Y(f,) und der Halbperiode a, 
herangezogen (#, ist jetzt Abschaltstelle von y), so gelten die Ungleichungen 0, =o,=o, 
also a,=Za,=a, und für #, erhalten wir die Schranken 


» 


die Strecke P,P, enthält den Nullpunkt im Innern. 


Für tz #, ist (1) durch das System (6) fortsetzbar, usw. In wiederholter Anwendung 
dieses Verfahrens erhalten wir für (1) eine spiralförmige Phasenkurve, welche bei >» von 


1) oder von innen (für ,<Y,<0) gegen die periodische 
Lösung des Systems I, konvergiert. Die Steuerungsfunktion Z(f) hat abzählbar unendlich 
viele einfache Nullstellen #,, t,, f, Aus (40), (41) und den weiteren analogen Ungleichungs- 
reihen folgt, daß die Differenz je zweier konsekutiver Nullstellen gegen @a(>0) strebt. Die 
Lösungen mit (35) und (39) erfüllen also (35) und die Trefftzsche Forderung des Abklingens 
der Störungen. 


außen (für 


Wir behaupten, daß Gl. (1) nur diese stabilen Lösungen besitzt. 


1. Für o=0 existiert keine periodische Lösung, der Fall a—=0 (unendlich große Frequenz) 
wurde ausgeschlossen. Die Bewegungen in I, sind instabil. Für Y,>0 ist (28) nicht er- 


füllt, es sei also - z*s Y,<0. Dann ist (1) durch (6) mit dem Anfangswert P, (X, =0, 


Y,<0) für £>0 fortsetzbar. X(f) nimmt zunächst von Null an monoton ab, erreicht für 
Y(#)=0 ein endliches Minimum in P,(X<0, Y=0); anschließend verläuft die Phasenkurve 
in der Halbebene Y>®, X(f) nimmt zu. Zur Zeit t=f,(>0) wird die Steuerungsgerade 
(hier die Y-Achse) in P,(X (£,)=0, Y(t,) >0) erstmalig erreicht. Es ist Y(t,)<!Y,|, da P, 


näher am Punkt Fe 0) liegt als P,. Über P, hinaus wird (1) durch (5) fortgesetzt, usw. Für 


{> strebt die Phasenkurve nach Null, bleibt also nicht stets im Innern von (35) mit r >0. 
2. Es sei o>®. Nach Il. gibt es in diesem Fall keine periodischen Lösungen. 


a) Y, mit »><Y,<0 sei die Anfangsgeschwindigkeit einer Bewegung, die für 


—R 
1+ 
t=0 im Phasenpunkt P,(X,,Y,) mit X,+0Y,=0 beginnen möge. Damit diese durch (1) 
beschrieben werden kann, muß nach (29) P, im Außengebiet von £,—=0 liegen. Da &,=0. 
eine Kurve zweiter Ordnung ist und im Ursprung die Y-Achse tangiert, schneiden alle Ge- 
raden @, die durch den Ursprung gehen und im Quadranten X, >0, Y,<0 der Nullebene 
liegen, diesen Kegelschnitt. Es sei also P, im Bereich £, >0 gewählt; (1) ist dann zunächst 
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durch (6) fortsetzbar bis f,, der ersten auf t=0 folgenden Schaltstelle, deren Existenz analog 
[Z, p- 323] sich beweisen läßt'*). Der Punkt P, (X (t,), Y (f,)) auf der Geraden Z—=0 liegt dem 
Nullpunkt näher als P,; P, und P, werden aber nicht notwendig durch diesen- getrennt. 
Liegt der Nullpunkt zwischen P, und P, und liegt überdies P, im (29) entsprechenden Be- 
reich für die Abschaltstelle von y, so ist (1) über #, hinaus durch (5) fortsetzbar; es gibt eine 
Schaltstelle tf,, usw. Die Möglichkeit der Fortsetzung ist jedoch höchstens endlich oft gegeben; 
die Folge der absoluten Beträge der Ordinaten Y(t,) dieser Schaltstellen hat nämlich keinen 
positiven Häufungspunkt (es müßte sonst eine periodische Lösung existieren, gegen welche 
die Bewegung durch P, konvergierte). Andererseits nimmt diese Folge monoton ab, es gibt 
daher ein erstes Y(t,) derart, daß der Punkt P,(X (t,), Y(t,)) auf der Steuerungsgeraden, 
jedoch nicht im Bereich (29) oder dem ihm entsprechenden einer Abschaltstelle liegt. Die Be- 
wegung ist hier wegen Y(f,)=F0 nicht beendet, sie hört aber auf, für t > #, der Differential- 
gleichung (1) zu genügen. 


b) Es sei Y, mit O<Y,< 


R 
Vi+e gegeben. Für #’< « ist eine reelle Bewegung für 


t>0 nicht durch (1) zu beschreiben. Falls #° > a ist (d.h. {,=0 eine Hyperbel), läßt sich (1) 
fortsetzen, wenn Y, derart gewählt wird, daß P,(X,,Y,) im Bereich {,<0 liegt. Es sei 
also P, innerer Punkt des Hyperbelastes Z{,=0 im zweiten Quadranten der Nullebene:; die 
durch P, gehende Steuerungsgerade schneidet diese Kurve. Analog a) tritt dann nach end- 
licher Zeit t, eine Schaltstelle P, auf, für welche die Bewegung zwar nicht beendet, aber 
dureh (1) nicht fortsetzbar wird. 


3. Es verbleiben die Lösungen mit o<0, Y, >0. 
Liegt P,(X,, Y,) im Innern von &,=0, ist also (29) 
erfüllt, dann ist (1) für alle £>0 definiert. Nach den 
Ausführnngen zu Bild 3 liegt P, in einem endlichen 
Bereich der Nullebene und (1) ist durch (6) zunächst 
für t>0 fortsetzbar. Wegen Y, >0 nimmt X(f) zu; 
für endliches £, >0 schneidet die Phasenkurve erst- 
malig die Steuerungsgerade in P, (X (t,) >0, Y(t,) >0). 
Die Phasenkurve liegt für O<t<t, z. B. ganz im 


Innern eines Kreises um (2, 0) mit einem Radius, der 


gleich dem Abstand dieses Punktes von P, ist. Über P, Bild 6. 

hinaus — einer Abschaltstelle von y, welche außerhalb 

&,=0 liegt — setzen wir (1) durch (5) fort. Der nächstfolgende Schaltpunkt P, (Bild 6) hat 
die Koordinaten X (t,)<0, Y(t,)<0'). Von P, an gilt unser früher bewiesener Satz, die 
Lösung konvergiert gegen die periodische Lösung in Io, gehört also von f, ab zur Klasse der 
bereits genannten stabilen Lösungen. 


Es ist noch anzumerken, daß in den zuletzt betrachteten Anfangswerten P, in 3. eine 
zweite stabile Lösung von (1) beginnen kann. Hierfür ist allerdings {=0 Abschaltstelle von y. 
Die zweite Phasenkurve durch P, besteht zunächst aus einem Bogen derjenigen Schar, welcher 
PP, in Bild 6 angehört; die weitere Fortsetzung ist dann evident. Für beide Lösungen stellt 
sich der gleiche Endzustand ein, nämlich die eindeutig vorhandene periodische Lösung in %%». 
Ebenso gibt es eine zweite Lösung durch die für eine Aufschaltstelle zulässigen Anfangs- 
werte P, im Fall (2b) mit #° > «u und o >0; aber auch diese Bewegung ist nur bis zu einem 
endlichen f fortsetzbar. Zur Vermeidung solcher Mehrdeutigkeiten ist daher erforderlich, 
das Gebiet £,<0O in (29) und das entsprechende bei einer Abschaltstelle als Anfangswerte 
auszuschließen. 378 


1) Vgl. G. D. Birkhoff: Dynamical Systems. Americ. Math. Soc. Publ. Bd. 9, New York 1927, 
S. 123 bis 128. 


“) Vgl. Birkhoff in loc. eit. '?). 
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Natürliche Schrittweite bei numerischer Integration 
von Differentialgleichungssystemen. 


Von L. Collatz in Karlsruhe. 


Die Brauchbarkeit der Interpolationswerfahren, die sich zur numerischen 
Integration gewöhnlicher Differentialgleichungen sehr bewährt haben, hängt 
stark von der Wahl der benutzten Schrittweite ab. Es gibt eine aus- 
gesprochene „günstigste Schrittweite“. Es wird eine Methode angegeben, 
nach der man bei einem beliebigen System gewöhnlicher Differentiul- 
gleichungen 1. und 2. Ordnung eine „natürliche Schrittweite“ bestimmen 
kann. Die bei der Rechnung tatsächlich benutzte Schrittweite sollte in 
ihrem Verlauf dem Verlauf der natürlichen Schrittweite möglichst an- 
gepaßt werden. 


1. Einleitung. 


Zur numerischen Integration gewöhnlicher Differentialgleichungen 1. und 2. Ordnung 
sind in neuerer Zeit die nach Adams und Störmer benannten Verfahren mehr und mehr 
angewandt worden; von diesen Verfahren zeichnen sich die Interpolationsverfahren gegenüber 
den Extrapolationsverfahren durch eine wesentlich höhere Genauigkeit bei nur mäßig ver- 
größerter Rechenarbeit aus. Es zeigt sich nun, daß die Durchführung der Interpolations- 
verfahren empfindlich von der Größe der zugrunde gelegten Schrittweite abhängt; es gibt bei 
vorgelegter Differentialgleichung oder bei einem System von Differentialgleichungen eine 
ausgesprochen günstigste Schrittweite A. Rechnet man nämlich mit kleinerer, etwa der halben 
Schrittweite, so verläuft die Rechnung zwar glatt und ohne Schwierigkeiten, aber man hat 
wegen der größeren Anzahl von Schritten zu dem gleichen Kurvenstück praktisch die doppelte 
Rechenarbeit zu leisten. Rechnet man aber mit der doppelter Schrittweite, so nehmen die 
höheren Differenzen nicht genügend stark ab, die Iterationen konvergieren zu langsam, man 
muß dreimal und unter Umständen noch öfters das Schema durchrechnen, ehe die Zahlen 
zum Stehen kommen, und die letzten mitgeführten Differenzen zeigen ein so starkes Schwanken, 
daß man zu ihrer Genauigkeit nieht mehr das nötige Vertrauen haben kann. 


In einem solchen Falle ist es stets besser, die Schrittweite zu verkleinern, etwa auf 
die halbe Sehrittweite überzugehen (der Übergang auf die halbe Schrittweite ist rechnerisch 
bequem ausführbar')). Der Bestimmung der in dieser Weise beschriebenen günstigsten Schritt- 
weite, die weiterhin als „natürliehe Schrittweite* bezeichnet wird, soll die folgende Arbeit 
dienen. Bei Kenntnis der natürlichen Schrittweite kann man dann die bei der Rechnung be- 
nützte Schrittweite diesem Verlauf jeweils durch Schrittverdoppelung bzw. Schritthalbierung 
anpassen. 


Ferner kann die Bestimmung der natürlichen Schrittweite dazu dienen, bei Integration 
eines Differentialgleichungssystems von vornherein die ungefähr benötigte Schrittanzahl und 
damit die erforderliche Rechenarbeit abzuschätzen °). 


Die natürliche Schrittweite bezieht sich auf das laufende Verfahren, dagegen nicht — 
und darauf sei ausdrücklich hingewiesen — auf die Berechnung eines Anfangsstückes, bei 
der die Iterationen in der Regel viel schlechter konvergieren. Man rechnet beim Anfangs- 
stück zweckmäßig mit einer zwei- bis dreimal so kleinen Schrittweite als beim laufenden 
Verfahren. 


Soweit die Beweise ähnlich denen der Arbeiten von G. Schulz’) verlaufen, konnten sie 
hier kurz gefaßt werden. Auch bei Verwendung anderer Integrationsverfahren, z. B. des 
Runge-Kutta-Verfahrens ist es nicht ratsam, stets mit gleicher Schrittweite h zu rechnen. 
Man kann die hier ermittelte „natürliche Schrittweite“ auch als Anhaltspunkt für die Wahl 
der Schrittweite bei diesen anderen Verfahren verwenden. 


1) L. Collatz-R. Zurmühl: Beiträge zu den Interpolationsverfahren der numerischen Integration von 
Differentialgleichungen 1. und 2. Ordnung, Z. angew. Math. Mech. Bd. 22 (1942), S. 42 bis 55. Dort sind die Verfahren 
zusammengestellt und verschiedene Fragen der praktischen Handhabung besprochen. 


2) Die Arbeit entstand im Anschluß an größere numerische Rechnungen, die am Institut f. prakt. Mathem. d. 
Techn. Hochschule Darmstadt ausgeführt wurden. Bei diesen Rechnungen wurde die hier vorgeschlagene Bestimmung 
der „natürlichen Schrittweite‘“ erprobt. 


3) G. Schulz: Z. angew. Math. Mech. Bd. 12 (1932), S. 44 bis 59 und Bd. 14 (1934), S. 224 bis 234. 
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2. Eine einzelne Differentialgleichung 1. Ordnung. 


In dieser Nummer sind einige bekannte Tatsachen zusammengestellt, auf Einzelheiten 
bei den Verfahren kann hier nicht eingegangen werden. 
Bei der Differentialgleichung 


genüge f(f, x) einer Lipschitzbedingung 


wobei in praktischen Fällen gewöhnlich X durch den größten Betrag der Ableitung J., !uner- 


halb des in Betracht kommenden Bereiches ersetzt wird: K=k,, mit k= 4 . Es seien 
bereits Näherungswerte x; für die exakten Werte «(t;) der Lösungsfunktion »(f) an den 
Stellen für öä=...n— 2, n—1,n berechnet. Bei fester Schrittweite ist ;=t,+ih. Der 


r . 0 . . 
Wert , wird durch Verbesserung eines Ausgangswertes , nach der Iterationsvorschrift 


gewonnen, wobei die üblichen Abkürzungen 


gebraucht sind. Dann gilt’) 
1. Die bei dem Verfahren nötigen Iterationen zur Ver- 
besserung der Werte von #„,;, konvergieren, wenn Kh<3 ist. 
2. Es läßt sich eine Fehlerabschätzung durchführen: der 
Fehler «,— x (t„) kann abgeschätzt werden durch 


15 


wobei ce, und e, von n unabhängig sind. Das Anwachsen des 
Fehlers kann also nicht stärker sein als das Anwachsen der + 
Potenz z". Die Größe z hängt nur von dem Produkt K'h ab 
(Bild 1). Es ist 


3 (2? —1) 
2+4z+1 | 


Für Kh>3 geht z>x. Die Zahl z darf natürlich nicht zu 
groß sein, z. B. schon z=1,3 (das entspricht Kh=0,%6) ist gila ı. Abhängigkeit der Größe 2 
reichlich groß. von Kh. 

Bei Verwendung eines anderen Verfahrens, z. B. des ge- 
wöhnlichen Adamschen Interpolationsverfahrens®), sind die Verhältnisse genau dieselben, 
nur tritt bei 1. an Stelle Kh < 3 die Bedingung Kh< 8/3, und bei 2. ist z aus einer anderen 
Kurve abzulesen ®). 

Für kleine Werte von h ist bei Benutzung der Formel (2.3) angenähert 


Z=1+Kh. 


Kh=W()= 


| 
| 


| 
- 
| 
| 


Einführung der „natürlichen Schrittweite“. Als wesentlich für die Konvergenz der Iterationen 
(d. i. gleichbedeutend mit der genügend starken Abnahme der Differenzen im Differenzen- 
schema) als auch für das Anwachsen der Fehlerschranke bei dem einzelnen Schritt erscheint 
also das Produkt Kh. Hat man nun bei einer Lösungskurve stark veränderliche Werte von k, 
so wird man, um längs der ganzen Kurve mit ungefähr denselben Verhältnissen zu arbeiten, 
die Schrittweite h nicht konstant lassen, sondern ebenfalls variieren, am bequemsten je nach 
Erfordernis zur doppelten bzw. zur halben Schrittweite übergehen, so daß kh längs der ganzen 
Kurve von derselben Größenordnung ist. So kann man den Wert 


konst. 


4) Der Gedankengang der Fehlerabschätzung ist bei G. Schulz3) durchgeführt. Für die Gl. (2.3) findet man 
das Ergebnis in Z. angew. Math. Mech. Bd. 22 (1942), S. 47. 

5) G. Sehulz: Formelsammlung zur praktischen Mathematik. Berlin und Leipzig 1937, Sammlung Göschen, 
Bd. 1110, S. 118. 

6) Die Kurve ist dargestellt Z. angew. Math. Mech. Bd. 22 (1942), S. 47. 
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als „natürliche Schrittweite* einführen. Der Wert von € kann 3 | | 7] 
nicht allgemein festgelegt werden; bei großen Genauigkeits- | | | 
ansprüchen wird man € kleiner wählen als bei Überschlags- a Ai 

rechnungen. Als Anhaltspunkt sei für C der Wertebereich I | | 


0,05 bis 0,15 genannt. x 


Beispiel: Die Differentialgleichung 


| 
| Ä 
| 
hat bei der Anfangsbedingung x (0) — 1 die Lösung | 
| I 
Hier ist k= Bild 2 zeigt die / 
Lösung x und die Größe e als Funktionen von f. Die Lösung rn | 
besitzt bei £=1 einen Pol, die Schrittweite muß bei An- ET 
näherung an die Stelle *=1 kleiner und kleiner gemacht mn 7 2, 3 4 
werden in Übereinstimmung mit dem Verlauf der natürlichen Bil 3; Beispiel für netirliche 


Schrittweite, die für >1 gegen den Wert O strebt. Schrittweite. 


3. Systeme von Differentialgleichungen 1. Ordnung. 
Bei dem System von m gewöhnlichen Differentialgleichungen 1. Ordnung 


tritt an Stelle der einen Lipschitzkonstanten K eine ganze Matrix 8 von m’ Lipschitz- 
konstanten 


m’ 


die definiert sind durch das Bestehen der Ungleichungen 


m 
Für praktische Zwecke kann man die K;, aus den Beträgen der partiellen Ableitungen be- 
rechnen, setzt also 


Konvergenz der Iterationen. Die Iterationen bei Berechnung der Näherungen «;,,;, an 
der Stelle t„;, werden wie bei (2.3) nach der Vorschrift ausgeführt 


v 


Es ist also 


m 
+1) (v) (v) (v—]) < (v—]) 


durch obere 


Hier kann man das Ungleichheitszeichen fortlassen, wenn man die «; 
Schranken X”) ersetzt. 
Faßt man die X’) „+1 als Komponenten eines Vektors x” auf, so gilt 


Der Grenzvektor 


€ 


. | 
| | | 
| 
| | | 
| 
€ 
4 
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charakteristische Zahlen hat, die dem Betrag nach sämtlich kleiner als 1 sind. Der größte 
Betrag der charakteristischen Zahlen von 8 sei K, dann lautet also eine hinreichende Be- 


dingung für die Konvergenz der Iterationen 
Kh<B3. 


existiert und die Konvergenz der Iterationen ist damit gesichert, wenn die Matrix 


Fehlerabschätzung. Es werden nun die Fehler eingeführt 
Dann gilt, wenn man von (3.4) unter Fortlassen der Indizes (v) und (r +1) die Gleichung 


(u-)+ fi (,xi.. At 
tn-ı 


abzieht und die Koeffizienten 


einführt: 


(3.2). 


= —1 
Auf der rechten Seite kann die erste eckige Klammer mit Hilfe der Lipschitzkonstanten ab- 
geschätzt werden, während die zweite den Quadraturfehler C;, , darstellt. Für den Quadratur- 
fehler C;,„ erhält man die Schranke (Fehlerglied der Simpsonregel)’) 


max 


wobei die vierte Abteilung für praktische Zwecke gewöhnlich durch den vierten Differenzen- 
quotienten ersetzt wird. Somit erhält man aus (3.5) 


m 


Schreibt man in (3.6) das Gleichheitszeichen, so erhält man für obere Schranken X;,„ der 
Größen &;,„, das System von Differenzengleichungen 


m 


2 

Ku, Be KrılXunsı-el+ C; . . . . . . (3.7). 
=] 


Hiermit ist im Prinzip die Fehlerabschätzung geleistet; kennt man die Fehler der Anfangs- 
werte, so kann man rekursiv obere Schranken für die weiteren Fehler berechnen. Man kann 
mit Hilfe der charakteristischen Zahlen der Matrix $£ eine Fehlerschranke in geschlossener 
Form aufstellen. Die GlIn. (3.7) stellen ein inhomogenes System von Differenzengleichungen 
dar. Eine spezielle Lösung des inhomogenen Systems ist X,,„—=konst.—=d;, wobei die d; aus 


m 
> 


zu berechnen sind. Die Determinante der K;; ist dabei als von Null verschieden voraus- 
zusetzen. Zu der so erhaltenen speziellen Lösung des inhomogenen Systems ist die allgemeine 
Lösung des homogenen Systems von (3.7) zu addieren. Die Lösung wird in üblicher Weise 
in der Form angesetzt 

X,n=Bi:2"; 
damit erhält man 


?),G. Kowalewski: Interpolation und genäherte Quadratur, Leipzig und Berlin 1932, S. 72, 


| 

> 


220 Collatz, Natürliche Schrittweite bei numerischer Integration Bd.» Nr. 4 Aug. 196 


Die Determinante dieses für die B; homogenen Gleichungssystems stellt gleich Null gesetzt 
eine Bestimmungsgleichung für z dar, deren absolut größte Wurzel z verwendet wird; denn 
es soll an Stelle der allgemeinen Lösung des Systems nur eine majorisierende Lösung auf- 
gestellt werden. Die Lösung ist dann nur brauchbar, wenn z>1 ist. Führt man die 


Funktion ein 


0 


so kann man (3.5) in der Gestalt schreiben 
m 


B;:W(e)=h K;ıBı. 
Dieses Gleichungssystem besagt: W(z) ist eine charakteristische Zahl der Matrix hs. Da 
W(z) für z2>1 monoton wächst und der größte Wert von z zu verwenden ist, kommt auch 
die größte charakteristische Zahl K der Matrix ft in Frage. 
Da W(e) für 1<z< nur Werte zwischen 0 und 3 annimmt, muß 


We)=hK<3 


sein. Das ist die gleiche Bedingung für die Schrittweite h, wie sie sich aus der Konvergenz- 
betrachtung für die Iterationen ergab. 

Zusammenfassend stellt man fest, daß bei einem System von Differentialgleichungen 
dieselben Verhältnisse herrschen wie bei einer einzelnen Differentialgleichung: man hat 
lediglich an Stelle der einen Lipschitzkonstanten K hier die größte charakteristische Zahl X 
der Matrik 8 der Lipschitzkonstanten zu setzen. Alles andere, sogar die Funktion W (z) 

konst. 
bleibt genau erhalten. Man wird daher auch hier entsprechend (2.6) die Größe h, = in 
als natürliche Schrittweite einführen, wobei k die größte charakteristische Zahl der Matrix 
mit den Elementen k;,, ist. Man nimmt also jetzt als Lipschitzkonstanten nicht mehr die 
Maximalwerte längs der ganzen Lösungskurve, sondern verwendet als K;, die Beträge der 
: A an der Stelle, an der man die natürliche Schrittweite bestimmen will und 

vergleicht das vorgegebene System von Differentialgleichungen mit einem System, bei dem 


Ableitungen 


_ konstante Werte haben, und zwar gerade die gleichen Beträge wie das 


vorgelegte System im herausgegriffenen Augenblick. Als Vergleichsmittel dient das An- 
wachsen der Fehlerschranke, also Hie Größe 2. Die hier gegebene Definiton von h, enthält 
natürlich noch eine gewisse Willkür; man könnte z. B., wie es sonst in der Theorie der 
Differentialgleichungen oft geschieht, an Stelle der Matrix (3.2) nur m Lipschitzkonstanten 


die Ableitungen 


* 


einzuführen. Dann tritt an Stelle der größten charakteristischen Zahl K der Matrix 8 ein- 
fach die Summe der K;. Die numerischen Rechnungen zeigen jedoch, daß man bierbei die 
Verhältnisse nur sehr viel gröber erfaßt als bei Benutzung der charakteristischen Zahlen. 


.4. Systeme gewöhnlicher Differentialgleichungen 2. Ordnung. 
Bei dem System 
mögen die gegebenen Funktionen f; den Lipschitzbedingungen genügen 
m 


Es werden Näherungswerte x, und x, für die exakten Werte x(t,) und &(t„) an der Stelle 
t„—=t,+nh berechnet nach dem Interpolationsverfahren der zentralen Differenzen (Berechnung 
der (»+ 1)-ten Näherung aus der v-ten nach): 


2 
| 
| 
m 
* K 7 “| 
= K; — x; | 
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(„+l) 
=2%,n - Ag n+1—o | 


(v+ 1) (9) 
+h 
Dabei sind die Konstanten 


1 10 


fi(tu+1-.) und fr Kaas bedeuten den exakten Wert von (4.1) an der Stelle —=#,;1_. und 
seinen im Verlauf der Rechnung dafür erhaltenen Näherungswert in der v-ten Näherung. 


Konvergenz der Iterationen. Entsprechend zu 3 erhält man hier für die Änderungen 
und 
ne 


+1) v—1) 


-(#+]1) 


h 


Sehätzt man die rechten Seiten mit Hilfe der Lipschitzbedingungen (4.2) ab, so folgt 


m 


—1) —]) 


Faßt man obere Schranken und der Größen und als Komponenten eines 2 m- 
dimensionalen Vektors v'”) auf 


so genügt v(») der Gleichung 
vr) — A 
wobei A die Matrix ist: 


h 
4 Lm m 4 Mmm 


Ist K, die dem Betrag nach größte charakteristische Zahl der Matrix A, welche nach einem 


Satz der Algebra, vgl. Ziffer 6, positiv ist, so ist also die Konvergenz der Iterationen ge- 
sichert, wenn 


ist. Sie ist nicht mehr gesichert, wenn Ky den Wert 7, unimmt, wenn also 


3 
det (a e) == 


ist (E ist die Einheitsmatrix). Eine kurze Umformung ergibt 


t 
1 
t m 
m 
I=1 
r 
1 
t 
r 
1 
a Tom Mm 
h 
.. 
| 
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3 Im 1 


3 
det (a e) 


Man kann diese Determinante auch auf die Gestalt bringen 
det E)—=0, 


wobei eine m-reihige quadratische Matrix mit den Elementen ist. 
Fehlerabschätzung. Es gilt für die exakten Werte 


wobei D;„ und D/„ abschätzbare und als bekannt anzusehende Quadraturfehler bedeuten. 
Zieht man diese von (4.3) ab (unter Fortlassen der Zeiger (v) und (r+1)), so erhält man für 
die Fehler 


bei Abschätzung der Differenz der Funktionswerte mit Hilfe der Lipschitzkonstanten: 


2 m 


(4.5). 


2 m 
—Ln-ıl sh (M;ı ı-e| + Liı | Din| | 
= 

Diese Zeilen bleiben richtig, wenn man die Absolutstriche fortläßt, Gleichheitszeichen an 
Stelle der Zeichen = setzt und an Stelle der &;,„, £;,n obere Schranken E;,„, Z;,n für sie 
schreibt. Dann entstehen Differenzengleichungen, die außer einer speziellen Lösung der in- 
homogenen Gln. (4.5) (als die man im allgemeinen E,,„=konst.=E;, Z;, „= konst. =Z,; ver- 
wenden kann), noch den zu (4.5) gehörigen homogenen Gleichungen genügende nach einem 

Potenzgesetz anwachsende Lösungen haben. Daher wird der Ansatz gemacht: 


Ei, n 4A; 2”, Z,, n = B; g" . . . . . . . . . . . (4.6). 


Es ergeben sich dann verschiedene Werte z, und für das Anwachsen der Fehlerschranke ent- 
scheidend ist das z mit dem größten Betrage. Beim Einsetzen von (4.6) in (4.5) erhält man 


2, m 

m 

i=1 
Mit den Abkürzungen 
= Wiz) 


wird daraus: 
m 


A; B;W(e) (A, + Bı 
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Diese 2m homogenen Gleichungen für die A; und B; haben nur eine nichttriviale Lösung, 
wenn ihre Determinante verschwindet. Sie lautet: 


Aus dieser Gleichung ist z zu berechnen. Die dem Betrage nach größte Wurzel z be- 
stimmt dann gemäß Gl. (4.6) das Wachstum der Fehlerschranken. Die tatsächliche Berechnung 
von 2 aus (4.8) ist im allgemeinen eine mühsame Aufgabe. Beschränkt man sich jedoch auf 
kleine Schrittweiten h, also h<1, so kann man mit dem Ansatz 

z=1+Kh 
in (4.7) und (4.8) eingehen; man erhält 


und (4.8) geht in eine gewöhnliche Säkulargleichung über, K ist dann die größte charakte- 
ristische Zahl der Matrix 8 mit nichtnegativen Elementen: 


Zur Abrundung soll noch auf einen Zusammenhang der Matrix 8 mit der Matrix 4 bei der 
Konvergenzuntersuchung der Iterationen hingewiesen werden. Die Fehlerabschätzung versagt, 
wenn z gegen © rückt. Für z>» geht nach (4.7) 


We)>7; 


Setzt man diese Werte in (4.8) ein, so unterscheidet sich die Determinante nur durch den 


m 
konstanten Faktor (3) von der Determinante (4.4). 


Natürliche Schrittweite. Bei veränderlichen Werten der partiellen Ableitungen Hi und 


Ai vergleicht man das vorgelegte System an einer herausgegriffenen Stelle t mit einem 


solchen, bei dem diese Ableitungen konstant sind und setzt: M;, = se, Lı= 

das Anwachsen der Fehlerschranke (4.6) (und zwar genau genommen für das Vergleichs- 
system) ist nach dem Vorangehenden das Produkt K'h entscheidend. Die größte charakte- 
ristische Zahl K der Matrix (4.9) ist dann eine Funktion der herausgegriffenen Stelle f. Es 
liegt daher nahe, als natürliche Schrittweite den Ausdruck 


konst. 
h=- K 


einzuführen, wobei über den Wert der Konstanten dasselbe gilt wie am Schluß von Ziffer 2. 
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5. Systeme von Differentialgleichungen erster und zweiter Ordnung. 

Liegt ein System von m Gleichungen zweiter und p Gleichungen erster Ordnung vor. 
. . O1), 
so setzt man y;—=&,„+; und hat dann ein System von m-+p Gleichungen 2. Ordnung, au! 
das sich die Theorie von Nr. 4 anwenden läßt. Schreibt man die zugehörige Determinante 
(4.8) auf, so stehen in den letzten p Spalten nur Nullen bis auf die Zahlen — O(z) in der 
Diagonalen; man kann daher die letzten p Spalten und Zeilen fortlassen und hat beim An- 

satz z=1-+Kh wieder K als größte charakteristische Zahl der Matrix 8 zu nehmen: 


m+n» 
m-+p 
Zeilen 


m 
Zeilen 


m Spalten p Spalten m Spalten 


oder in leichtverständlicher symbolischer Abkürzung 
| 
| 
E 0 0 


Diese Matrix enthält alle bisherigen als Spezialfälle. 


6. Berechnung der charakteristischen Zahl vom größten Betrag. 


Eine m-reihige quadratische Matrix 


(6.1) 


hat m charakteristische Zahlen A,, K,,... AK; das sind diejenigen Werte von X, für die das 
homogene Gleichungssystem 
m 
Naıc=Ke; 

(oder in Matrizenschreibweise Ar—=Kr) eine von identisch Null verschiedene Lösung 
X=(&,,%,...%m) hat. Die bis auf einen konstanten Faktor bestimmte Lösung heißt der zu 
der betreffenden Zahl gehörige Eigenvektor. Bei uns hat nun die Matrix (6.1) nur positive 
(oder verschwindende) Elemente. Für solche Matrizen gilt der Satz: Diejenige von den m 
charakteristischen Zahlen, die den größten Absolutbetrag hat, ist positiv und heißt die 
Maximalwurzel von X (sie kann auch eine mehrfache Nullstelle von det (A — KE)=0 sein). 
Es gibt zu ihr einen Eigenvektor, dessen Komponenten nichtnegativ und nicht alle Null sind 
(Satz von Frobenius)?). 


Zur praktischen Bestimmung von K wähle man einen beliebigen Ausgangsvektor u mit 
nur positiven Komponenten #,,%,,... %, und bilde den Vektor v—A u mit den Komponenten 


m 


v= Ur. 


3) G. Frobenius: Sitzungsber. Preuß. Akad. d. Wissensch. Math.-phys. Klasse, Berlin 1912, 1. Halbband, S. 456 
bis 477, insbes. S. 457. 
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. 
ie Quotienten «;—=— aus den einzelnen Komponenten von d und den entsprechenden von u 


ı 
werden nicht für alle © den gleichen Wert haben; man sieht aber leicht, wo ınan durch Ab- 
ändern der einzelnen Komponenten erreichen kann, daß die Größt- und Kleinstwerte der 
Quotienten «; enger zusammenrücken. Es gilt der Satz, daß die Maximalwurzel AK stets 
zwischen dem Größtwert und Kleinstwert der Quotienten liegt ?). 


Min u; <= K= Max u;. 


Zudem kommt es nicht auf sehr hohe Genauigkeit än, es genügt für die natürliche Schritt- 
weite, daß man die erste geltende Stelle sicher hat. Bei etwas Übung läßt sich die Be- 
stimmung der größten charakteristischen Zahl auf eine geltende Stelle sehr rasch durchführen. 


Beispiel: Vorgelegt sei die Matrix 
09 21 
2 0 


Weiß man gar nichts über den Eigenvektor, so wird man etwa von u=(l,1,1) ausgehen; 


dann wird 
b=Au=(25; 7,2; 2,6), 


bei den neuen Komponenten ist die mittlere viel größer als die anderen beiden, also wird 
jetzt etwa u—=(1,2,1) gewählt, das ergibt v=Au= (3,8; 8,1; 4,6). Die Quotienten sind 3,8; 
4,05; 4,6. 

Zur Herabminderung des letzten Quotienten 4,6 wird man die letzte Zahl 1 bei u er- 
höhen und zur Erhöhung des ersten Quotienten 3,8 wird man die erste Zahl 1 bei u etwas 
verkleinern. So erhält man u= (0,96; 2; 1,12); v = (3,912; 8,184; 4,576). Quotienten: 4,075; 
4.092; 4,086; das genügt für den vorliegenden Zweck vollauf; es ist also K=4,08, wobei die 
Dezimale 8 natürlich unsicher ist. Mit Sicherheit weiß man nur 


4,075 < K<4,0R. 


7. Rechnungsgang. 


Vorgelegt sei ein System (3.1), (4.1) oder (5.1) von gewöhnlichen Differentialgleichungen 


weit diese gebraucht werden), entweder durch Differenzieren oder durch Bildung von Differenzen- 


Sso- 


il. und 2. Ordnung. Zunächst bestimmt man die partiellen Ableitungen „A (bzw. 


quotienten (z. B. bei nur empirisch gegebenen Zusammenhängen). Die Werte dieser partiellen 
Ableitungen werden für eine Anzahl Zeitpunkte t=1", {”... ermittelt, nämlich für die Werte 
von f, für die man die natürliche Schrittweite bestimmen will. Mit den Beträgen der Ab- ° 
leitungen kann man die Matrix $ aufstellen, (3.2) bzw. (4.9) oder (5.2), indem man 


of: i M; ofi | 


setzt. Für jeden Zeitpunkt t=1) wird die größte charakteristische Zahl K der zugehörigen 
Matrix $? nach Nr. 6 ermittelt. Dann gibt u _i den zeitlichen Verlauf der natürlichen 
Schrittweite. Der Wert von C hängt von der gewünschten Genauigkeit ab und ist bei größeren 
Genauigkeitsansprüchen kleiner zu nehmen. Als grober Anhaltspunkt sei der Wert C=0,1 
genannt. 


Die Schrittweite, mit der man tatsächlich rechnet, wird man der natürlichen Schritt- 
weite möglichst anpassen, wobei man im ersten Teile der Rechnung eher mit kleineren 
Schrittweiten arbeitet und am Schlusse der Rechnung größere Schritte nehmen darf, weil 
die dabei entstehenden Ungenauigkeiten sich nicht mehr so stark auswirken können wie an- 
fangs entstandene Fehler. Die hier bestimmte natürliche Schrittweite bezieht sich nicht auf 
die Anfangsiteration für die Ausgangswerte, die stets besser mit kleinerer Schrittweite vor- 
genommen wird. 373 


®) L. Collatz: Mathem. Zeitschr. (1942) erscheint demnächst. 
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Über gleichwahrscheinliche Aufteilungen. 
Von H. Hadwiger in Bern. 


Ein statistisches Schema, das im Gebiet der angewandten Mathematik ver- 
schiedenartige Realisierungen erfährt, wird durch die gleichwahrscheinliche 
Zerlegung einer positiven Größe in positire Teile geliefert. In dieser Arbeit 
wird die Wahrscheinlichkeit, daß die Teilgrößen alle eine gegebene Schranke 
nicht überschreiten, einer eingehenden Studie unterzogen und ferner der 
Erwartungswert für die größte Teilgröße ermittelt. 


1. Problemstellung und Überblick. 


Die Summe von » unabhängigen natürlichen Zahlen sei 


Für die Summanden x, sei jede der Zahlen 0, 1,2,...s gleichwahrscheinlich, oder anders aus- 
gedrückt, jede Lösung der diophantischen Gl. (1) sei mit der gleichen Wahrscheinlichkeit zu 
erwarten. Wir wollen in einem solchen Fall von „gleichwahrscheinlichen Auftei- 
lungen“ sprechen. Es muß mit besonderem Nachdruck auf diese Voraussetzung hingewiesen 
werden, um einer Verwechslung mit den in der matlıematischen Statistik häufig untersuchten 
„statistischen Aufteilungen“ vorzubeugen. Diese Aufteilungen werden so erzeugt, 
daß s Elemente auf n Plätze verteilt, und die einzelnen x, in (1) als Besetzungszalilen der 
Plätze interpretiert werden. In diesem Falle sind, wie bekannt, die verschiedenen Auftei- 
lungen (1) keineswegs mehr gleichwahrscheinlich, sondern müssen mit den Gewichten 


1\s s! 


versehen werden, vorausgesetzt, daß die Zuteilungswahrscheinlichkeiten für die einzelnen 
Plätze alle gleich sind. 
Wir stellen nun die folgenden Aufgaben an die Spitze unserer Untersuchungen: 
A) Wie groß ist die Wahrscheinlichkeit X,„(&), daß rt ausfällt für alle 


B) Wie groß ist der Erwartungswert ©, für =, wo x, die größte der Zalılen 
(r=1,2,...n) bezeichnet ? 
Während die Behandlung dieser Fragestellungen bei statistischen Aufteilungen auf erhebliche 
Schwierigkeiten zu führen scheint'), zeigt es sich, daß bei gleichwahrscheinlichen Aufteilungen 
eine befriedigende Lösung erzielt werden kann. 


Wir werden das oben gestellte Problem in kontinuierlicher Gestalt formulieren, und dann 
die exakten Lösungen entwickeln. 


Wir stellen zunächst kurz die Hauptergebnisse zusammen: 
Für die gesuchte Wahrscheinlichkeit ergibt sich 


wobei als Hilfsfunktion 


1 n 


gewählt wurde. Ferner erhält man für den gesuchten Erwartungswert 


0 


also eine bemerkenswert einfache Formel. Nach einem bekannten Theorem gilt für große n 
die asymptotische Formel 


n=„(C+logn) 


wo C die Eulersche Transzendente € — 0,577215 ... bezeichnet. 


ı) Vgl. z.B. Titeica, Serban: Sur un problöme de probabilitcs. Bull. Math. Phys. Eeole polyteehn. Bucarest. 
Bd. 10 (1939), S. 57 bis 64. 
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Die Hilfsfunktionen (4) stellen symmetrische Verteilungen dar, die in einem gewissen 
Fragenkreis der Wahrscheinlichkeitsreehnung und Fehlertheorie eine wichtige Rolle spielen. 
Nach einer Formel von A. Sommerfeld?) stellt 


das (n — 1) dimensionale Volumen des Schnittes des » dimensionalen Einheitwürfels mit einer 
senkrecht zu einer Würfeldiagonalen stehenden „Ebene“ im Abstand « vom Würfelmittel- 
punkt dar. 
Nach einer bekannten Methode kann diese Volumformel mit Hilfe des Diskontinuitäts- 
faktors von Dirichlet 


cos et dt— 


- 


hergeleitet werden’). 
Nach G. Poölya*) gilt die Grenzwertformel 


sn 


eine Konvergenz, die in jedem endlichen Intervall gleichmäßig stattfindet. 


Die Funktionen (4) können durch Polynome dargestellt werden, und zwar ist nach einer 
Formel von Laplace°) 


D„(x) = 


(n 


wo beidemal das gleiche Vorzeichen gewählt werden soll®). 


Durch geeignete Interpretation des beim gestellten Problem vorliegenden Sachverhaltes 
kann die Lösung (3) auf Grund der Formel (7) von Sommerfeld gewonnen werden. Wir 
schlagen jedoch hier einen andern Weg ein, um auf rein analytischem Wege zum Ziel zu 
gelangen. 


2. Integraldarstellungen. 
Wenn wir in (1) #,=£,s setzen, so gilt 


Wir fassen nun die £, als kontinuierliche Variable auf, und geben damit dem oben gestellten 
Problem die neue endgültige Form, in welcher die gesuchte Funktion als geometrische Wahr- 
scheinlichkeit, der gesuchte Erwartungswert als geometrischer Mittelwert neu definiert wird. 

Wie üblich stützt sich die Wahrscheinlichkeitsdefinition auf ein Maß, das durch ein 
Riemannsches Integral gegeben wird. Mit Rücksicht auf die wesentliche Voraussetzung, 
wonach jeder möglichen Aufteilung (11) das nämliche Gewicht zukommen soll, haben wir 
als Maß für eine kontinuierliche Menge von Aufteilungen das Integral 


zu setzen, wo sich die Integration über die genannte Menge zu erstrecken hat. Die Existenz 
des Integrals wird dabei vorausgesetzt, was einer entsprechenden Beschränkung der zur 
Messung zugelassenen Aufteilungsmengen gleichkommt. 

Bezeichnet J,„ das durch (12) erklärte Maß sämtlicher Aufteilungen (11), so definieren 
wir als Wahrscheinlichkeit dafür, daß eine Aufteilung zu einer vorgegebenen Teilmenge M 
gehört, den Integralquotienten 


2) A. Sommerfeld: Eine besonders anschauliche Ableitung des Gaußschen Fehlergesetzes. Boltzmann Fest- 
schrift. Leipzig 194, S. 848. Vgl. weiter: V. Brun: Norsk. Mat. Tidskr. Bd. 14 (1932), S. 81 bis 92; P. Hertz: Math. 
Ann. Bd. 74 (1913), S. 153 bis 203. 

°*) J. M. David: Zur Dirichletschen Methode des Diskontinuitäts-Faktors. Diss. Zürich 1918, S. 45 bis 49. 

4) G. Pölya: Bereehnung eines bestimmten integrals. Math. Ann. Bd. 74 (1913), S. 204 bis 212. 

5) Theorie analytique des probabilites, S. 169. 

%, In der von G. Pölya (vgl. das Zitat in Fußnote 4) angegebenen Formel fehlt infolge Druckfehlers der 
Faktor 
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7, 

Das Maß aller Aufteilungen 

das wir mit J„(£) bezeichnen wollen, ist leicht zu finden. Nach (12) gilt zunächst 


wo die Integration über das Gebiet 


zu erstrecken ist. Indem man zuerst &, ,—» festhält, und erst nach erfolgter Integration 
über die übrigen Veränderlichen noch über n integriert, gewinnt man die Rekursion 


woraus sich die Formel 
. . . . . . . . . . . . . (17) 


ergibt. Insbesondere ist also 


1 


Für die Wahrscheinlichkeit, daß für eine Aufteilung (11) die Bedingung 


erfüllt ist, erhält man nach (13) das Integral 
3 


wo die Wahrscheinlichkeitsdichte 9, (£) durch 
J 
n 


gegeben ist. 
Für die Wahrscheinlichkeit, daß für eine Aufteilung (11) die Bedingung 


erfüllt ist, erhält man nach (13) das Integral 


das über das Gebiet 
1-8 —,—...-&-,SE 


zu erstrecken ist. 
Wir werden aber zur Berechnung dieser Wahrscheinlichkeit nicht von dieser Integral- 


darstellung ausgehen, sondern diese als Lösung einer Funktionalgleichung entwickeln. 
Bezeichnet £&, die größte der Zahlen &,,&,,...£,„, so ist die Wahrscheinlichkeit dafür, 


daß für eine Aufteilung (11) die Bedingung 
erfüllt ist, offenbar durch 


gegeben. Nach geläufiger Schlußweise erhält man für den Erwartungswert ©, von &, 
1 
9,=1-—| . . . . . (22). 
0 


3. Eine Funktionalgleichung. 


Wir werden nun für die Wahrscheinlichkeiten X,(£) eine Funktionalrekursion herleiten. 
Es sei n=2. Wir betrachten zunächst die Wahrscheinlichkeit dafür, daß für eine Auf- 
teilung (11) die Nebenbedingungen 
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MAIN 
‚nr 


erfüllt sind. Die Wahrscheinlichkeiten für die einzelnen Bedingungen sind: 
a) (nach (19)), 
(man dividiere in a) durch £), 


Nach dem Produktsatz ist die oben bezeichnete Wahrscheinlichkeit 
(an — X,._, 


Wenn man nun die Bedingung a) fallen läßt, so ergibt sich die Wahrscheinliehkeit X,,(£) 


nach dem Summensatz 
ı 


— 1)\ 
0 
oder mit anderer Berücksichtigung der Hilfsfunktion » (£) 


A, (£ = (n d . . . . . . . . . (23). 


Die Funktionalgleichung (23) stellt eine Integralrekursion dar, die es ermöglicht, sämt- 
liche Funktionen durch fortgesetzte Integrationen zu gewinnen. Selbstverständlich muß die 
Ausgangsfunktion gegeben sein. Offenbar ist 


0(-»<i<I) 


(24). 


Auf diese Weise erhält man 


X, 


» usw. Durch eine geeignete Transformation der Funktion X,,(£) kann die Funktionalrelation (23), 
die noch eine ziemlich undurchsichtige Gestalt aufweist, erheblich vereinfacht werden. Wir 
machen den Ansatz: 


Die Relation (23) wird dadurch transformiert in 
1 
2 
n)= \ Bu —n+ 1)a 


1-5 


a) 
<: 
e) 
14), | 
15), 
16), | 
IS). 
( 
(- 
und 
1 ( 1S!E<o) 
-2+65-3# ( 
0 (-»<:=5) 
1) 
r | 
ıf- 
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Durch die weitere gleichzeitige Transformation von Haupt- und Integrationsvariablen 


gewinnt man 


x+l 


So erhalten wir für die Hilfsfunktionen ®, eine bemerkenswert einfache Integralrekursion, 
die übrigens eine anschauliche Interpretation zuläßt: 
Der Funktionswert ®,(x) ist der Integralmittelwert von ®,„_,(x), der über die sym- 
metrische Nachbarschaft von von —1 bis erstreckt wird. 


Wie man aus (24) in Verbindung mit (27) folgert, ist die Ausgangsfunktion für die Re- 
kursion (28) 


1 
— 

0 


Mit Hilfe der Funktion (29) kann die Funktionalrelation (28) noch auf die folgende Form ge- 
bracht werden 


4. Eine Darstellung der Lösung. 


Wir wollen nun für die sich aus der Integralrekursion ergebenden Funktionen ®, (x) 
eine explizite Darstellung finden, um so mit (27) zu einer geschlossenen Formel für die ge- 


suchte Wahrscheinlichkeit X,„(£) zu gelangen. Es bezeichne g„(t) die Fouriersche Trans- 
formierte 


welche sicher existiert, da ®,„(x) außerhalb des endlichen Intervalles |«|<n verschwindet, 
wie man zurückgreifend auf die Darstellung (27) leicht verifiziert. 

Durch Transformation der Funktionalrelation (30) erhält man nach dem bekannten 
Faltungssatz ’) die algebraische Beziehung 


Wie man leicht nachrechnet, ist aber 
sin t 

so daß aus (32) folgt 
sin t\r 

Nach der Inversionsformel für die Fouriersche Transformation ergibt sich nun 


Mit (27) und (35) ist das erste Hauptresultat (3) und (4) erreicht. 
Wie man aus der Integraldarstellung (35) abliest, gilt 


Aus dieser Symmetrie kann noch eine merkwürdige Funktionalgleichung abgeleitet werden. 
Wie man auf Grund von (27) leicht verifiziert, ergibt sich als Folge von (36) 


?) Vgl. etwa G. Doetseh: Theorie und Anwendung der Laplace-Transformation. Berlin 1937, S. 162. 
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Da für das Intervall 


das Argument 


ausfällt, liefert die Funktionalgleichung für dieses Intervall die Funktionswerte, nämlich 


1 


5. Eine Rekursion. 

Wir wollen nun für die durch (22) definierten Erwartungswerte ©, eine Rekursion her- 
leiten, aus welcher sich dann der explizite Ausdruck (5), also das zweite Hauptresultat mühe- 
los ergibt. 

Im folgenden sei stets n=3. Wir setzen 


wo also nach (22 
1 


1 
0 
oder unter Verwendung der Darstellung (27) 

1 

12 
0 0 


ist. Durch die Transformation t=Fa ergibt sich 


1 


2 
h=z | | (in cos 


00 
oder wenn wir die Hilfsfunktionen 
= | (sin eos 2a —na)da 
0 
(a) = \ (sin 
0 
einführen, die Darstellung 


0 
Naheliegende Abschätzungen des Integrals (43) sind 


wo die erste für große a, die zweite dagegen für kleine « in Anspruch genommen werden kann. 
Durch zweimalige partielle Integration von (45) gewinnt man unter Beachtung von (46) 


(46), 


(47). 


Durch Anwendung der trigonometrischen Additionstheoreme ergeben sich für die Hilfs- 
funktionen (43) und (44) die Rekursionen 


und 


1 1 
=! 
on, 
te- 
9) 
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x) 
ze 
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en 
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Durch Verbindung der erhaltenen Relationen (48) und (49) gewinnt man jetzt 


sin a sin (n — 2)a 1 


Andererseits ergibt zweimalige Ableitung von en 
—2)a) (sina"'"— (ld) . 
Setzen wir den in (50) geschriebenen Ausdruck in (51) ein, so erhalten wir 


wo zur Abkürzung 
= In cos (n - cos(n —3)a+ (4 sin «a sin (n — 2) al (sin . (53), 
oder nach einfacher trigonometrischer Umformung 


1 1 1 
gesetzt wurde. Verwenden wir den für die zweite Ableitung gewonnenen Ausdruck (52) in 
der Formel (47), und drücken die vom Teil (54) herrührenden bestimmten Integrale durch die 
Hilfsfunktionen (4) aus, so resultiert 


2 1 | An- 2 .. 
Nun ist nach der Formel (10) von Laplace 
1 
so daß wir die Rekursion 
1 1 1 
. . ö . . (57) 
erreichen. Nach (39) ist aber 
zu setzen. So ergibt sich aus (57) die gesuchte Rekursion 


Da die Rekursion um zwei Stellen zurückgreift, müssen zwei Ausgangswerte gegeben sein. 
Im Hinblick auf die Voraussetzung n =3 sind das offenbar 


die sich aus (22) unter Verwendung von (24) und (25) berechnen lassen. 
Es gelingt nun, die sich auf Grund der Rekursion (58) ergebende Folge in expliziter 
Gestalt anzugeben. Wie man mühelos direkt verifizieren kann, ist nämlich 


Damit haben wir das zweite Hauptresultat (5) gewonnen. 


6. Tabelle der Erwartungswerte. 
Wir geben nun eine Zusammenstellung der ersten 20 Erwartungswerte 9,: 


n n ©, 
1,000 000 6 0,408 333 11 0,274534 16 0,211 296 
| 0,750 000 70,370 408 12 | 0,258 601 17 | 0,202 327 
0,611 111 0,339 732 13 0,244 636 18 0,194 173 
0,520 833 9 0314330 14 0.232354 19 0,186 723 
0,456 667 10 0,292 897 20 | 0,179 887 363 
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KLEINE MITTEILUNGEN 


Ergänzungen zu: Kamke, Diiferentialgleichuagen: 
Lösungsmethoden und Lösungen. 


Herr Görtler veröffentlicht nachfolgend eine 
Reihe von Differentialgleichungen mit ihren Lösun- 
gen als Ergänzung zu 'dem kürzlich erschienenen 
Buch von E. Kamke, Differentialgleichungen: 
Lösungsmethoden und Lösungen, Leipzig 1942. Da 
die im dritten Teil dieses Buches gegebene Samm- 
lung von Lösungen von Differentialgleichungen 
ihre Aufgabe um so besser erfüllen wird, je voll- 
ständiger sie ist, hat Herr Görtler angeregt, 
in der Zamm eine Rubrik für die Veröffentlichung 
solcher Ergänzungen zu dem Buch einzurichten. 
Wir halten diese Anregung für sehr glücklich und 
bitten daher alle Leser der Zeitschrift, durch Ein- 
sendung etwaiger von ihnen gelöster und noch 
nicht  veröffentlichter Differentialgleichungen zur 
Vervollständigung der Sammlung beizutragen. 
Jeder einzelne einmal bearbeitete Fall kann 
anderen von Nutzen sein und Arbeit ersparen. 

Einsendungen bitten wir an den Herausgeber der 
Zamm zu richten. Die Veröffentlichungen werden 
unter dem obigen Sammeltitel erfolgen. 


Tübingen/Dresden. Kamke. Willers. 


Ergänzungen. In der kurzen Zeit seit dem ersten Er- 
scheinen der Kamkeschen Sammlung von Differen- 
tialeleichungen nebst Lösungen (Vorausdruck 1940 
im Auftrage der Zentrale für wissenschaftliches Be- 
richtswesen der Luftfahrtforschung des Generalluft- 
zeugmeisters) hat sich der überaus große Nutzen 
einer solchen Sammlung für die Praxis erwiesen. 
Der Aufforderung, diese Sammlung möglichst bald 
um alle erreichbaren, noch nicht aufgenommenen 
Gleichungen mit Lösungen (vor allem aber auch 
um jene, die, der Allgemeinheit nichts nutzend, in 
den Schreibtischen unveröffentlicht schlummern) zu 
ereänzen, komme ich zunächst, um einen Anfang 
zu machen, mit einigen anspruchlosen Vertretern 
aus der Klasse der linearen Differentialgleichungen 
zweiter Ordnung nach. Ich habe die teilweise zu- 
fällig, teilweise im Rahmen bestimmter Frage- 
stellungen der angewandten Mathematik gzefunde- 
nen Ergebnisse nur so weit verallgemeinert, als 
mir dies aus Zeitgründen gegenwärtie möglich 
war. Immerhin umfassen viele der Beispiele Son- 
derfälle, welche bereits in der Sammlung aufge- 
nommen sind. Ich habe diese, insofern sie mir auf- 
gefallen sind, jeweils unter Anführung der be- 
treffenden Ziffer des Kamkeschen Buches !) an- 
gegeben. Die Anordnung der Gleichungstypen ent- 
spricht der Kamkeschen Systematik. 


+e?+(2ac+b)cosr 


+(2be—a)sinr 


a? b? 
cos2r—+ 


y=expier+asinz— bcosr} 
[C,+C;,\exp{— 
(Hillsche Differentialgleichung). 

[C, +6, t— 2(ar-+be*)y dr]. 


!)E.Kameke: Differentialgleiehungen: Lösungsmethoden 
und Lösungen. I. Gewöhnliche Differentialgleichungen. 
Leipzig 1942, Teil C. 


.....6) 
a 


Yy=exp = e*!?) mitv= 


wo Z, die allgemeine Zylinderfunktion der 
Ordnung » bezeichnet. 


. (4). 


1.Ist =a(a—a), , =b(1+a—.a), so ist 
y=expi—ar—bet} 
2.It ,=ala—a), b,=ab, so ist 
y"+ta+ber}y ta (4a); 
It ,=ala— a), 3, 
so ist 
y=exp{—ar+Be*} 
[C,+C,\expt@a— 
((4) ist hierin enthalten. Sonderfall:Kamke 2.34.) 
y=expfia\fdr} 
+C,\expt— 
(Sonderfall: Kamke 2.29 für f(x) = —g(r).) 
zy" y' 
+t{a+b,2 y=0 


1. Ist a =b(a —1)+B(a—2), 
c=c(b-+B), so ist 


„2 


+0, [20 exp] + “las 


(T). 


(Sonderfall: Kamke 2.136.) 


2. Ist a =b(a—1)+PBla-2), 
so ist 


J cx? 
Cı+C, | x“ — s-tdzl. 
(Sonderfall: Kamke 2.121.) 


3.Istt m =-aß, =—-B(b+P), so ist 


(Sonderfall: Kamke 2.141.) 


4.It „= —aß, 
c,=c(b-+B), so ist 


y=exp 


[a 


(Sonderfall: Kamke 2.122. 
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zy"+2(1+ br) y' 


| 


[C,+C,\exptr— 2ce*}da]. 
y | 
(9; 


+ z{a+bx)y 
+- fa(a—a— 1)+(alb 


3) + (10); 


y=zterdxic, Vader]. 
(Sonderfälle: Kamke 2.176; 2.179: 2.187; 2.196; 
2.198; 2.199; 2.202; 2.277; 2.281). 
ey" 
SCHEN. 
+ta, da? 


1. Für -+Bla-au), 
lautet die Lösung: 


Yy fx 


2. Für 
d=c(b—P), 
lautet die Lösung: 


y=ı*exp - 
(Sonderfälle: Kamke 2.209, ferner natürlich (7)). 


a) 
+(a Na) | 


(12): 


— 


a2 y" +2 (f(@) + 2a) y’ | 

+ ta(a— 1) 2? +(a+ba)F(a)y 


(13); 


y=exp{—if!de) 
de). 
383 


Göttingen. H. Görtler. 


Berichtigungen zu: Kamke, Differentialglei- 
chungen, Lösungsmethoden und Lösungen. Seite 160, 
Zeile 9, statt (x) lies: fs(y). Seite 431, 2: 126a, 
statt des zweiten y” lies: u”. 


Kleine Mitteilungen Bd. 22 Nr.4 Aug. 1942 
Über die Lösung von Gileichungs- 

systemen durch Iteration. Im Folgenden 

sell eine hinreichende Bedingung für die Lösbar- 

keit eines Gleichungssystem mit mehreren Unbe- 

kannten durch (spezielle) Iteration mitgeteilt 

werden. Der Einfachheit halber wird sie nur für 


zwei Unbekannte formuliert werden, woraus ilr 
Wortlaut für mehr Unbekannte sofort ersichtlich 
ist. Das System habe die Form 


z=g(T,Y) 


Die durch hinreichend gute Ausgangsnäherungs- 


werte (z'®,y'®) für eine Lösung (a,b) und durelı 
die Rechenvorschrift 


bestinmten Wertefolgen y'”’) konvergieren be- 
kanntlich gegen die Lösung (a, b), falls 


Dabei sind My. My,Nx,Ny obere Schranken von 
|dy 
bzw. 
dr °y|’ |dz 
Bereich der z-y-Ebene, der sämtliche Näherungen 
y'’”) und die Lösung (a,b) enthält). 
Diese Bedingung erweist sich aber vielfach als 
sehr ungünstig. Bei Gelegenheit einer Fehlerab- 


in einem konvexen 


schätzung bei der numerischen Integration von 
Differentialgleichungen (Z. angew. Math. Mech. 


Bd. 12 [1932], S. 44—59) war vom Verf. bereits eine 
günstigere angegeben worden (a. a. O., S. 56-57). 
die aber wegen der speziellen, auf das dort be- 
handelte Problem zugeschnittenen Formulierung 
offenbar wenig bekannt geworden ist. 

Aus (2) folgt nach dem Mittelwertsatz 


rv 
Op 
(a 
(4), 
( Y 


wobei sämtliche partiellen Ableitungen an einer 
Zwischenstelle der geradlinigen Verbindungsstrecke 
der Punkte (2 y”’ =») und (z°”, y'”) zu nehmen 
sind. Daraus folgt weiter 


_ 


M,„ My N. Ny die obigen Bedeutungen 
Besitzt nun die Matrix nichtnegativer Ele- 


wobei 
haben. 
mente 


My) 
\ N x N Y 


nur charakteristische Zahlen, das heißt Wurzein 
der Gleichung 
uw M 
Nx Ny—n 
die dem Betrage nach kleiner sind als 1, dann 
konvergieren die Reihen 


1) Vgl. Fr. A. Willers: Methoden der praktischen 
Analysis, Berlin und Leipzig 1928, S. 167. — H. v.Sanden: 
Praktische Analysis, 2. Aufl., Leipzig und Berlin 1923, 
S. 160 bis 162. 


5). 
( 
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und erst recht die Reihen 


L (er — + (2 — lim .r 


(4), 
= | 
und zwar gegen die Lösung (a, b). 

Zum Beweis wird der folgende Hilfssatz ver- 
wendet, der bereits a. a. 0. S. 52-53 formuliert 
wurde: Bestehen zwischen den Gliedern von k 
unendlichen Reihen 


mit positiven Gliedern für alle vr >N die k linearen 
Ungleichungen 


und besitzt deren Matrix E = {c,.} nur nichtnega- 
tive Elemente sowie charakteristische Zahlen, die 
dem Betrage nach kleiner als 1 sind, dann konver- 
gieren die k Reihen sämtlich. 

Für die wichtigen Fälle von zwei oder drei Un- 
bekannten ist die Probe, ob die Bedingung erfüllt 
ist. ohne große Mühe durchzuführen, zumal die 
charakteristischen Zahlen nur sehr roh bestimmt 
zu werden brauchen. Bei ihrer Berechnung macht 
man bisweilen vorteilhaft davon Gebrauch, daß 
nach einem Satze von G. Frobenius die cha- 
rakteristische Gleichung einer Matrix mit positiven 
Elementen eine Wurzel hat, die reell, positiv. ein- 
fach und absolut größer ist als jede andere Wurzel 
(die sogenannte Maximalwurzel. Kommen ver- 
schwindende Elemente in der Matrix vor, so bedarf 
dieser Satz einer Abänderung: indessen kann ja 
immer erreicht werden, daß sämtliche obere Schran- 
ken positiv sind. 

Es mag noch erwähnt werden, daß auch 


N+N,<V<i. 8 


eine hinreichende Bedingung ist. 
schließt man von (5) aus weiter, daß 


gm < (Mx + M,) 
= (N. + N,) 
max — 27V]. 


Zum Beweis 


) 
(9) 


yr yr-» ) 


gilt und daher auch 


Daraus folgt dann die absolute Konvergenz der 
Reihen (7). 

Schließlich ist es auch für die Konvergenz der 
Iterationen hinreichend, wenn die Quadratsumme 
aller Elemente in der Matrix der Schranken für 
die Absolutbeträge der partiellen ersten Ableitun- 
gen .kleiner als 1 ist, wenn also bei zwei Unbe- 
kannten 


(10). 


ztMy + (M) 


ist 2). Ist nämlich die Summe der Quadrate aller 
Elemente einer (reellen) Matrix kleiner als 1, so 
folgt nach zweimaliger Anwendung der Schwarz- 
schen Ungleichung, daß die Folge der Potenzen 
der Matrix gegen die Nullmatrix konvergiert. Dann 
und nur dann sind aber die charakteristischen 
Zahlen sämtlich dem Betrage nach kleiner als 1. 
und man kann wie oben weiterschließen. 


Berlin. Günther Schulz. 384 


?) Für Systeme linearer Gleiehungen findet sich eine 
entsprechende Bedingung bei R. v. Mises und H. Polla- 
ezek-Geiringer: Z. angew Math. Mech., Bd. 9 (1929), 
S. 58 bis 77 u. S. 152 bis 164, und zwar auf S. 67. 


Zur Interpolation bei ungleichen Tafel- 
abständen. Die Aufgabe, einen Näherungswert 
der Funktion f(x) an einer Stelle x zu berechnen, 
falls für die Argumente .... 7, die zuge- 
hörigen Funktionswerte Yo Yır Y„ zegeben 
sind, wird am einfachsten mit der Lagrangeschen 
oder der Newtonschen Interpolationsformel gelöst. 
Vielfach zeigt sich nach Ausführung der Inter- 
polation, daß man noch weitere Funktionswerte 
in Rechnung zu ziehen hat. Während in einem 
solchen Falle bei Verwendung der Lagrangeschen 
Interpolationsformel die ganze Rechnung von 
vorne durchzuführen ist. kann man bei der New- 
tonschen Interpolationsformel die bisherige Rech- 
nung benutzen. weshalb sie für die Praxis die 
größere Bedeutung hat. 

Neben dieser Forderung, bei Erweiterung der 
Interpolation die bisherige Rechenarbeit verwen- 
den zu können, die man an eine praktische Inter- 
polationsformel stellen muß, ist es außerdem er- 
wünscht, die Änderung des bisherigen Näherungs- 
wertes, die durch Hinzunahme eines weiteren 
Funktionswertes erzielt wird. überschlagsmäßig 
rasch abschätzen zu können und außerdem diese 
Änderung möglichst schnell zu erhalten, vor allem 
falls sehr viele Interpolationen durchzuführen sind. 

Bei Verwendung der allgemeinen Newtonschen 
Interpolationsformel 


y=y+ 2](z — 
— 2) (e—2)+:-- 

— 


pflest man die darin auftretenden „Steigungen“ 
die symmetrische Funktionen von 
7, sind z. B. aus 
1 Yv Yv +1 


Ip — 


zu berechnen !). 

Diese Interpolationsformel hat. auch wenn der 
Näherungswert an einer x, benachbarten Stelle be- 
rechnet werden soll und die Zahlen #4. -... 
in zweckmäßiger Weise auf beiden Seiten von x 
liegen, vor allem bei großen Abständen der Argu- 
mente 2», den Nachteil. daß die Steieungen 
im allgemeinen mit wachsendem 
sehr klein, dagegen die Produkte (r — r,) (r — 
_,) entsprechend groß werden, daß 
man ferner schlecht übersehen kann, in welcher 
Weise die Hinzunahme eines neuen Funktions- 
wertes das bisherige Ergebnis der Interpolation 
beeinflußt. 

Sind die Abstände der x, alle gleich. +, 


h, und setzt man ° =f, so erhält man (1) 


h 
unter Benutzung absteigender Differenzen in der 
für die praktische Rechnung bequemeren Gestalt 


0 1 


.9 


tt +2) 
1-2-3 


tt +2)---(t+n--1) 


+ 


wobei 4,” das 1. Glied der vten Differenzenfolge 
von YoYı---%Yn bedeutet, A,” = — 4," 


1) Vgl. z. B. F.A. Willers: Methoden der praktischen 
Analysis, Berlin und Leipzig 1928, S. 58 ff, 
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Liegt x in der Nähe von 2, 9 — hK<e<a+!h, 
so sind die Koeffizienten der Differenzen in (3) 
absolut genommen alle <1 und der Betrag des 
bei Erhöhung des Grades der Interpolation hinzu- 
tretenden Gliedes ist absolut genommen jedenfalls 
kleiner als der Absolutbetrag der zugehörigen 
Differenz. 

Es sei hier für den Fall ungleicher Abstände 
der Argumente eine andere Form der Interpola- 
tionsformel angegeben, die sich an die Form für 
gleiche Abstände anschließt. Sie eignet sich vor 
allem, falls in ein und derselben Tafel sehr viele 
Interpolationen durchgeführt werden müssen, z. B. 
bei Einargumentinterpolationen in einer Tafel von 
mehreren Veränderlichen. Bei Verwendung ent- 
sprechender Hilfstafeln erlaubt sie den durch Hin- 
auftre- 


zunahme eines weiteren Funktionswertes 
tenden Beitrag schnell abzuschätzen und zu be- 


rechnen. 
Wir bilden 


2m -1) 


P == 1 
Iy— 


und damit die verallgemeinerten „Differenzen“ 


D,® 


Dy+ 


u=1,2,3,.-- 


also das folgende Schema 


und erhalten aus (1) 


2) 
+Ds 2) (2. — 2) 
+D,® (X —2)(2 (e—2,) (6). 
En-ı) 
D,® 
Es gilt nämlich 
D,® 0% 


(£y—Ly +1) Ey Lynn) 


Für »—1 ist das offensichtlich: 
— 
Ip — Tv tr 


dann durch 


(7) folgt 


Nach (2) ist 


 ; 


vollständige Induktion. 


(Ey 


— 
D,® 
1 
— 
1 
D,®+% 


Setzen wir 


dy = 


1 1 


‚=23-- 


wobei 


1 
da, =4,_ı pP. 


so erhalten wir aus (6) 


y=4+ Di’ t+Do® 


In P,“”"und in a, kommt x bzw. ? nicht mehr 
vor, so daß man sich, falls sehr viele Interpola- 
tionen auszuführen sind, eine Hilfstafel für P," 
und die Folge «a, anfertigt und damit mit Hilfe 
des Rechenschiebers oder einer geeigneten Rechen- 
maschine ?2) schnell den neuen Beitrag ausrechnen 
kann. 

Sind alle Tafelabstände gleich, so sind alle 
P,® =1, also = 4,", a,=» und (10) geht 
unmittelbar in (3) über. 

Bei Interpolation aus einer Tafel wird man als 
... Argumente wählen, die abwechselnd 
auf beiden Seiten der Interpolationsstelle z ge- 
legen sind. Wir schreiben für solche Fälle das 
Differenzenschema in folgender Weise: 


2) vgl. F. Emde: Fortlaufende Rechnungen auf der 
a 2. Instrumentenkde. Bd. 56 (1936), S. 185 
is 188 
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wobei also 


em 


mit Pevr-ı = 1 
Iy-u— 2 5 


x 


[6777] 
=> 
2 
a2 em 
29 +1 v +1 
? 2 


= Yr — 
2 


am 
-Dar+i . 
2 


D, =D;v-ı —P 


v 
2 


Liegt dann z zwischen z, und r,. so wird man 
die Anzahl und Reihenfolge der in der Umgebung 
des Punktes z gelegenen Tafelwerte z, je nach 
ihrer Dichte und Lage in der einen oder anderen 
Weise für die Interpolation wählen. Nimmt man 
so erhält man aus (1) 


+ [2 2_,)(2— 2) (2 — 
+92 (2 — 2) z,) 


(= — 


Wegen der Symmetrie der Steigungen ist z. B. 


0) 


Dı 
2 


- 


analog bei den anderen Steigungen. 

Demnach erhalten wir bei dieser Einteilung unter 
Heranziehung der 2rn + 1ten Differenzen als 
Näherungswert 


Yo + 


—% 
Lo — 


(?_,—2%) (2_1-2ı) 


(r — (r 
-1 


+ (x = (r 
(r (2 
— T_n4e) 

— Lu) 

— ı) (r n” Tu) 

7 


eninD 
Dı \ 
(Fan 
— 2 — 
2) (x n — 
Hier sind also die in einer Ziekzacklinie um r ge- 


legenen Differenzen benutzt. Man wird / 


wieder als neue Veränderliche einführen. Ebenso 
kann man weitere Formeln bilden je nach der 
Lage der Tafelwerte bezüglich x, z. B. solche 
analog zu den Formeln von Gauß, Stirling und 
Bessel für gleiche Tafelabstände. Dabei treten 
hier an Stelle der einfachen Mittel 


1 an +n 1 au) 


die allgemeineren Mittel 


ı +P, Dı 
2 


bzw. 


au) 


). 


Für die Interpolation in manchen Tafeln mit un- 
gleichen Tafelabständen wird man mehrere solcher 
Formeln benutzei müssen, je nach der Lage und 
der Dichte der Tafelwerte in der Umgebung von vr. 

Diese Formeln erweisen sich vor allem bei der 
Interpolation einer Funktion von mehreren Ver- 
änderlichen durch Interpolationen nach je einer 
Veränderlichen als zweckmäßig, falls die Tafel für 
eine oder mehrere Veränderliche ungleiche Tafel- 
abstände besitzt, aber für alle tabulierten Werte 
der Veränderlichen die Einteilung nach einer Ver- 
änderlichen stets dieselbe ist. In diesem Falle 
stellt man sich für jede Veränderliche das Schema 
der P,„“ und bei Interpolation nach (10) oder einer 
der anderen Formen die Folge der a, bzw. der 
entsprechenden Größen auf. 

Für die in vielen Fällen ausreichende quadrati- 


und den meist in den Tafeln stehenden ersten 
Differenzen den Näherungswert 

1 

2 
er 3 
+1 


sche Interpolation liefert (10) mit 


(13) den Näherungswert 


1 (15). 
+1) 
(Liegt die Interpolationsstelle im Intervall 2,2, 
näher an z, als an z_,, so wird man im allge- 
meinen (14), sonst (15) benutzen.) Für die Skala 
der Tafelargumente berechnet man die Tafel der 
P. ferner für jedes Intervall die beiden Koeffizienten 

1 
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und kann den quadratischen Beitrag unter bBe- 
nutzung der meist in den Tafeln angegebenen 
ersten Differenzen z. B. mit einer Rechenmaschine 
mit Komplementschieber und Mehrfachmultiplika- 
tion rasch berechnen. 

München. 376 


J. Heinhold. 


über die Berechnung von beliebig lan- 
gen Aggregaten und „Faktorgregaten‘“ 
mit Hilfe eines einzigen einfachen Nomo- 
gramms. 1. Die Aggregate wollen wir annehmen 
in der Gestalt 

während wir unter einem Faktorgregat einen Aus- 
(lruck verstehen wollen von der Gestalt 


EM) 


so daß also der Logarithmus eines Faktorgregates 
eleich wird mit dem Aggregat der Logarithmen 
der Glieder. 

Dann hat Verf. bereits früher *) darauf hinge- 
wiesen, daß man (im Falle von lauter Pluszeichen) 
(lie Berechnung der zugehörigen Summen und Pro- 
dukte bei beliebiger Länge vornehmen kann dureh 
geeignete mehrmalige Verwendung eines einzigen 
Nomogramms für Addition. bzw. Logarithmen- 
Addition. Inzwischen hat Verf. bemerkt, daß man 
auch beliebig viele Glieder mit negativem Vor- 
zeichen zulassen kann, daß also auch die Berech- 
rung beliebiger Aggregate bzw. Faktorgregate auf 
diesem Wege ebenso leicht möglich ist. 

2. Nun zur Sache. Wir dürfen uns dabei auf 
die Berechnung von Aggregaten beschränken, da 
ja — wie schon erwähnt — der Logarithmus eines 
Faktorgregates gleich dem Aggregat der Glieder- 
logarithmen ist. 

Man nehme also eine gewöhnliche Netztafel für 
Addition und zeichne darauf links von und parallel 
zur y-Achse (Nullpunkte gleich hoch mit dem Ko- 
ordinatenanfang) je eine Skala für /(n) und 9 (@). 
sowie unterhalb von und parallel zur .-Achse 
(Nullpunkte senkrecht unter dem Koortdinaten- 
anfang) weiter Skalen für 

dabei sind die Skalenrichtungen gleich oder ent- 
gegen der Richtung der .„-Achse zu wählen, je 
nachdem das Zeichen + oder — gilt. und es 
können die letzteren Skalen ganz oder teilweise 
zusammengeleg‘ werden, wenn dies möglich ist. 
Die Skalen seien gegeben durch die Formein 

yllın), 

,=+19,(), 
Die Skaleneinheitsstrecke ist also überall dieselbe. 

Nunmehr gestaltet sich die Verwendung der so 
erhaltenen Tafel folgendermaßen: Man geht vom 
Punkt (0, 4.) horizontal und vom Punkt (.,.0) ver- 
tikal zum Punkt (,. 40). von da schräg (in Rich- 
tung const.) zum Punkt (0.40 + 
dann geht man allgemein vom Punkt (0. vo 
+...+:»_,) horizontal und vom Punkt (.r,,0) 
vertikal bis zum Punkt + %ı+...+7,_1) 
sowie von da schräg zum Punkt (0,w +7, +... 
+7,»_,+:%,). So kommt man zuletzt zum Punkt 
0,v)=(0,Wv+tXı+...+7,) und hat damit den 
Wert berechnet. 

3. Wenn es sich lediglich darum handelt. Ag- 
eregate von nur positiven Größen zu berechnen. 
so kann man — um Papier zu sparen — auch 
folgendermaßen vorgehen: Man superponiert eine 
eewöhnliche Additionstafel mit einer (ebenso ge- 
wöhnliehen) Subtraktionstafel, d. h. man zeichnet 
auf nicht zu engem quadratischem Netzpapier zwei 
nicht zu dichte Scharen von Geraden 


1) AWF-Mitteilungen 21. Jahrg. (1939), Heft 5, S. 77. 


3 €” 43 
Bill. 


Auf dieser Doppeltafel zeichnet man zwei Skalen 
für /(n) und 9,(5) wie vorhin, sowie bis zu n 
Skalen für 9,(E), ..., 9, (E®), diesmal 
sämtlich in gleicher Richtung mit der .-Achse., 

Die Verwendung dieses Rechenpapiers gestaltet 
sich dann ähnlicn wie soeben, mit dem einzigen 
Unterschied, daß die . hinter dem Komma ein 
Vorzeichen # bekommen, und man vom Punkt 
27,-,) zum Punkte W,% + x, 
+ 7,_, # dadurch gelangt. daß man ent- 
lang der ersten oder zweiten Geradenschar geht, 
je nachdem für das Glied .r, das Zeichen + bzw. 

- eilt. 

4. Zum Schlusse weisen wir darauf hin, daß man 
die hier zur Rede stehenden Funktionszusammen- 
hänge, besonders die relativ häufig auftretenden 
und darum wichtigen Faktorgregate, bisher immer 
in viel komplizierterer Weise dargestellt hat. Vgl. 
neben den Beispielen in der unter !) zitierten Note 


noch das Buch von H. Körwien,. Graphisches 
Rechnen (Nomographie), Deutsche Werkmeister- 
Rücherei, etwa 1940, Fig. 89—91, wo die Formel 
sehr umständlich mittels dreigliedriger (zusammen- 
gesetzter) Tafeln dargestellt ist. . 
7 
[7 
70 
Yı 
Bild 2. 


Münster i. W. Ludwig Neder. 38 
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Dr.-Ing. WILHELM CAUER, B. Apl. Prof. a. d. 
Teehn. Hochschule Berlin u. Leiter des Labora- 
toriums der Mix & Genest AG. Theorie der 
linearen Wechselstromschaltungen, 
1. Bd. XI +614 S. m. 426 Abb. Leipzig 1941, 
Akademische Verlagsgesellschaft Becker & Erler 
Kom.-Ges. Preis geb. 45 M. 

Für die Theorie und die Berechnung elektrischer 
Schaltungen waren die Kirchhoffschen Ge- 
setze und die komplexe Rechnung während vieler 
Jahre die einzige und auch ausreichende Grunid- 
lage. : Erst die Erfindung der elektrischen Sieb- 
ketten und die damit zusammenhängenden Pro- 
bleme der elektrischen Nachrichtentechnik stellten 
die Schaltungslehre vor neue Aufgaben, deren 
Lösung weitergehende Hilfsmittel erfordert. Wäh- 
rend man z. B. den Frequenzgang der Eigen- 
schaften einer gegebenen Schaltung mit den alt- 
bekannten Methoden leicht berechnen kann, reichen 
diese zur Lösung der umgekehrten Aufgabe nicht 
hin. Den ersten Schritt zur Bestimmung von 
Schaltungen, deren Betriebseigenschaften in ge- 
wünschter Weise von der Frequenz abhängen, tat 
Cauer auf der Grundlage eines von R.M. Foster 
1924 gezebenen Reaktanztheorems in der 1927 er- 
schienenen Arbeit über „die Verwirklichung von 
Wechselstromwiderständen vorgeschriebener Fre- 
quenzabhängigkeit“. Im Anschluß hieran hat 
Cauer in tiefgehenden Untersuchungen die all- 
gemeinen mathematischen Zusammenhänge zwi- 
schen den Schaltungen und. ihren Betriebseigen- 
schaften erforscht und für die praktische Berech- 
nung ausgewertet. Über einen Teil dieser Arbeiten 
und auf denselben Gegenstand gerichtete Unter- 
suchungen von 0. Brune, W. Brandt, 
H. Piloty und anderen berichtet das vorliegende 
Werk. 

Seinen Inhalt kennzeichnen die Kapitelüber- 
schriften: I. Aufgabenstellung und Beispiele. 
II. Berechnung der Eigenschaften gegebener 
Schaltungen. III. Vierpole. IV. Positive Funk- 
tionen und positive Matrizen. V. Reaktanz- 
theoreme. VI/VH. Wellenparametertheorie der 
Tiefpaß-Reaktanzfilter bzw. allgemeiner Filter. 
VII. Reaktanzvierpole mit vorgegebenen Betriebs- 
eigenschaften. IX. Frequenzweichen. X. AÄqui- 
valenz von Reaktanzschaltungen. 

Der Verfasser behandelt in knapper Form und 
in rühmenswerter mathematischer Strenge einen 
ausgedehnten Komplex von Problemen aus der 
Theorie der elektrischen Netze. Er beschränkt 
sich dabei keineswegs auf die allgemeinen theo- 
retischen Grundlagen, sondern dringt bis zur prak- 
tischen Berechnung von technischen Einzelauf- 
gaben vor, die durch zahlreiche Beispiele und An- 
leitungen erläutert wird. Allerdings ist das Buch 
nicht leicht lesbar. Der Verfasser sagt im Vor- 
wort, daß „an mathematischen Vorkenntnissen 
nicht mehr vorausgesetzt wird, als was man bei 
einem Durchschnittsabsolventen einer Technischen 
Hochschule voraussetzen darf“. Nach meiner 
Kenntnis der Verhältnisse dürften nur wenige aka- 
demisch vorgebildete Ingenieure über das Maß an 
Vorkenntnissen verfügen, die das Studium dieses 
Puches tatsächlich erfordert. Störend wirken die 
Unterbrechungen der Gedankenfolge durch einge- 
streute mathematische Exkursionen, die nur in 
losem Zusammenhang mit dem Thema _ stehen. 
Wer sich ungeachtet dieser Schwierigkeiten die 
Mühe nicht verdrießen läßt, das Buch gründlich 
durchzuarbeiten, wird reichen Gewinn davontragen. 

Die Berechnung der Filter hat der Verfasser 
etwas einseitig auf die Approximation nach 
Tschebyscheff abgestellt: es wäre er- 
wünscht, daß er daneben auch noch mehr auf 
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andere Arbeiten einginge. Ferner möchte man 
sich für die Neuauflage die Beseitigung von 
sprachlichen Amerikanismen, wie „Laplace-Ent- 
wicklung“, „Brune-Prozeß“, „Tschebyscheff-Ver- 
halten“ (statt „Laplacesche Entwicklung“ usw.) 
wünschen. Diese kritischen Bemerkungen sollen 
den hohen Wert des Cauerschen Brches nicht 
herabsetzen. 


Berlin. K. W. Wagener. 33 


Dr.-Ing. MAX SCHILHANSL, Die mittra- 
gende Breite bei der Kreisplatte mit 
radialen Rippen. {VDI-Forsehungsheft 411, 
Beilage zu „Forschung auf dem Gebiete des In- 
genieurwesens“, Ausgabe B, Bd. 12, Nov./Dez. 1941.) 
25 S. m. 13 Bildern und 9 Tafeln. Berlin 1941. 
VDI-Verlag G.m.b.H. Preis brosch. 5 M. 

Für die belastete Kreisplatte mit radialen Rippen 
werden Spannungen und Verschiebungen ermittelt. 
indem Kreisplatte und Rippen getrennt als Schei- 
benprobleme behandelt werden und in den Ver- 
bindungslinien gleiche Radialverschiebungen un:l 
Schubspannungen verlangt werden. Die Einleitung 
gibt ein vorzügliches Bild, wie das Problem ange- 
faßt worden ist. Weiter folgen umfangreiche für 
‚len angewandten Mathematiker beachtenswerte 
Berechnungen und Umformungen. Zum Schlusse 
wird ein der Rippenlage entsprechendes Sternsystem 
von T-förmigen Trägern angegeben mit gleichen 
elastischen Randeigenschaften: die Flanschbreite 
der T-Träger wird als mittragende Breite bezeichnet. 


Dresden. C. Weber. 39 


Dr. ©. SCHAEFER, Prof. a. dl. Universität Bres- 
lau, u. Dr. L. BERGMANN, Prof. a. d. Technischen 
Hochschule Breslau Grundaufgaben des 
physikalischen Praktikums, bearbeitet 
von Dr. W. Kliefoth, Lehrbeauftragter a. d. 
Universität und Technischen Hochschule Breslau, 
VI+218 S. m. 182 Abb. Leipzig und Berlin 1942, 
Verlag B. G. Teubner. Preis geb. 5.60 M. 

Das Erscheinen der vorliegenden Anleitung zum 
Physikalischen Praktikum ist von Lehrenden und 
Lernenden sicher mit größter Genugtuung begrüßt 
worden, da es beim praktischen physikalischen 
Anfängerunterricht eine allgemein empfundene 
Lücke ausfüllt. Angesichts des derzeitigen großen 
Mangels an Lehrkräften ist für einen erfolgreichen 
Unterricht im Physikalischen Praktikum eine gute 
Anleitung die unerläßliche Voraussetzung. Wir 
haben zwar für die Physiker und im besonderen 
auch für die Mediziner in den Büchern von West- 
phal bzw. Gerthsen-Pollermann sehr nützliche Ein- 
führungen ins Physikalische Praktikum. Doch 
fehlte es bisher an einer geeigneten Anweisung für 
alle anderen Gruppen, wie Chemiker und In- 
genieure., 

Das Buch bringt nicht nur eine genaue An- 
weisung zur Durchführung der Versuche, sondern 
auch zu jeder Aufgabe eine Einführung in ihre 
physikalischen Grundlagen. So bildet das Buch 
einmal eine sehr wertvolle Ergänzung zum eigent- 
lichen Lehrbuch und vermag außerdem die Unter- 
richtsassistenten außerordentlich zu entlasten. 

Die Aufgabenzusammenstellung entspricht etwa 
dem an den meisten Physikalischen Instituten 
üblichen Praktikum. Über die Auswahl im ein- 
zelnen kann man selbstverständlich verschiedener 
Auffassung sein. Doch gehören einige Aufgaben 
aus der Strömungslehre, sowie Aufgaben zur all- 
gemeinen Schwingungslehre, z. B. über erzwungene 
Schwingungen, den Übergang zur Resonanz, ze- 
koppelte Schwingungen, auch zum Aufgabenplan 
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für Ingenieure. Ihre Berücksichtigung in der 
nächsten Auflage würde daher den Wert des vor- 
trefflichen Büchleins nur steigern. 


Dresden. Stuart 8381 


Dr. phil., Dr.-Ing. e. h., Dr. phys. e. h. R. W. 
POHL, o. ö. Prof. a. d. Univ. Göttingen, Ein- 
führungin die Mechanik, Akustik und 
Wärmelehre. 3. u. 4 Auf. X +32 S. m. 
527 Abb. Berlin 1941. Verlag J. Springer. Preis 
geb. 15,80 M. 


Die beiden ersten Auflagen dieses ausgezeich- 
neten Lehrbuchs waren bereits in dieser Zeitschrift 
(Bd. 10 [1930] S. 419 und Bd. 12 [1932] S. 124) be- 
sprochen worden. Die günstige Beurteilung, die 
das Buch damals erfuhr, kann in verstärktem Mau 
für die neue Auflage gelten. Die Knappheit und 
Schärfe des Ausdrucks, die glückliche Wahl sinn- 
fölliger und einprägsamer Beispiele, die große An- 
schaulichkeit der Abbildungen, namentlich der für 
die Pohlschen Bücher charakteristischen Schatten- 
risse, die Vollständigkeit und Fülle des behandelten 
Stoffs sichern dem Buch auch weiterhin seinen 
Rang unter den besten Physikbüchern. Zu dem in 
den früheren Auflagen behandelten Stoff, der teil- 
weise umgearbeitet, teilweise, namentlich durch 
ein Kapitel über die Eigenschaften der festen Kör- 
per, erweitert worden ist, tritt in der neuen Auf- 
lage die Wärmelehre. Den Hauptinhalt bildet, nach 
einer kurzen Einführung der Grundbegriffe, die Be- 
handlung «der beiden Hauptsätze, wobei wiederum 
die Erzielung klaren Verständnisses für die physi- 
kalische Bedeutung und die technische Anwenil- 
barkeit der Sätze durch einfache Beispiele und an- 
schauliche Darstellung vollkommen erreicht wird. 
Die Einführung des Entropiebegriffs erfolgt im 
Anschluß an die Besprechuag nichtumkehrbarer 
Vorgänge als Maß der Nichtumkehrbarkeit, der 
zweite Hauptsatz wird im Anschluß an die Berech- 
nung des erößtmöglichen Wirkungsgrades einer 
Wärmekraftmaschine formuliert. Hervorzuheben 
wäre vielleicht noch als charakteristisch für die 
Darstellung die Verwendung des für die Technik 
wichtigen Enthalpie-Entropie-Diagramms. Neuartig 
in dem Buch ist — einem an der Phövsikalisch- 
Technischen Reichsanstalt eingeführten Gebrauch 
folgend — die Benutzung des Wortes „Kiloponil“ 
an Stelle von „kg-Kraft“, während das Wort „Kilo- 
gramm" die „kg-Masse“ bezeichnet. 


Berlin. E. Mosch. 389 


Ferner sind bei der 
Bücher eingegangen (ausführliche 
bleibt vorbehalten): 


Schriftleitung  folgen.e 
Besprechung 


Dr.-Ing. WILLIBALD LICHTENHELDT, Ein- 
fache Konstruktionsverfahren zur 
Ermittlung der Abmessungen von 
Kurbelgetrieben (VDI-Forschungsheft 408, 
Beilage zu „Forschung auf dem Gebiete des In- 
genieurwesens“, Ausgabe B, Bd. 12, Mai/Juni 1941). 
22 S. m. 71 Bildern. Berlin 1941, VDI-Verlag G.m. 
b.H. Preis brosch. 5 M. 


Dr. phil. nat. habil. MAX STECK, Dozent der 
Math. a. d. Techn. Hochschule München, Das 
Hauptproblem der Mathematik. XVI 
+ 232 S. Berlin 1942, Dr. Georg Lüttke Verlag. 
Preis geb. 10,60 M. 


Dr.-Ing. RICHARD EMIL KUTTERER, Balli- 
stik, mit einem Vorwort von Geh. Reg.-Rat Dr. 
C. Cranz. (Die Wissenschaft, Einzeldarstellungen 
aus der Naturwissenschaft u. der Technik, Bd. 97. 
herausgegeben von Prof. Dr. Wilhelm Westphal.) 
VII 208 Braunschweig 1942, Verlag Friedr. 
Vieweg & Sohn. Preis geb. 13,50 M. 


WERNER HEISENBERG, Wandlungen in 
den Grundlagen der Naturwissen- 
schaft. 6 Vorträge. 3. erw. Aufl. 95 S. Leip- 
zig 1942, Verlag S. Hirzel. Preis kart. 3,50 M. 


Dr. FRANZ BAUR, Prof. a. d. Universität Frank- 
furt/M.. Einführung in die Großwetter- 
forschung. (Math.-Phys. Bibl. Reihe I: heraus- 
zegeben von W. Lietzmann und W. Kliefoth, 
Bd. 88.) 2. verb. Aufl. 53 S. m. 14 Abb. im Text. 
Leipzig-Berlin 1942, Verlag B. G. Teubner. Preis 
kart. 1,20 M. 


Prof. Dr. PAUL RIEBESELL, Mathematik 
des täglichen Lebens. 52 8. Hamburg 
1942, Hanseatische Verlagsanstalt. Preis kart. 2M. 


THEODOR VAHLEN, Ballistik. 2. Aufl. neu- 
bearb. u. herauszeg, unter Mitwirkung von Alfred 
Klose. X + 267 S.m. 65 Abb. Berlin 1942, Verlag 
Walter de Gruyter & Co. Preis geb. 18 M. 


Prof. G. MAHLER, Physikalische Auf- 
geabensammlung. Neu bearb. von Prof. K. 
Mahler, Studienrat a. d. Schubartoberschule Aalen 
in Württemberg. 6. ver. Aufl. (Sammlung Göschen 
Bd. 243). 128 S. Berlin 1942, Verlag Walter de 
Gruyter & Co. Preis geb. 1,62 M. 


Dr.-Ing. P. WERKMEISTER, em. o. Prof. a. d. 
Techn. Hochschule Dresden. Vermessungs- 
kunde, Il. Stückmessung und Nivel- 
lieren. 7. Aufl. (Sammlung Göschen Bd. 468). 
165 S. m. 145 Fig. Berlin 1942, Verlag Walter de 
Gruyter & Co. Preis geb. 1.62 M. 


NACHRICHTEN 


Wissenschaftlicher Wettbewerb der Gesellschaft 
für Zeitmeßkunde und Uhrentechnik e.V. 

Die Gesellschaft für Zeitmeßkunde und Uhren- 
technik schreibt wie in anderen Jahren auch für 
das Jahr 1942/43 einen Wettbewerb für wissen- 
schaftliche Arbeiten aus. Zugelassen sind wissen- 
schaftliche Arbeiten von Wert aus den Gebieten 
der Zeitmeßkunde und Uhrentechnik. Ein festes 
Thema ist nicht vorgeschrieben. Für Preise steht 
ein Betrag von 2000,— RM zur Verfügung. Diesen 
Wettbewerb betreffende Anfragen sind zu richten 
an die Gesellschaft für Zeitmeßkunde und Uhren- 
technik e. V., Berlin SW 68, Neuenburger Str. 8. 
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Berichtigung zu der Besprechung des Buches 
von Dr. H. Blasius: Wärmelehre, Physikalische 
Grundlagen vom technischen Standpunkt. Z. angew. 


Math. Mech. Bd.22 (1942), S.62. „Herr Blasius 
weist mich darauf hin, daß er die Enthalpie i mit 
„Energiestrom“, nicht, wie es in meiner Besprechung 
infolge eines Schreibfehlers heißt, mit ..Wärmestrom“ 
bezeichnet. Meine Bedenken beziehen sich auf die 
Möglichkeit der Verwechslung dieser Bezeichnung 
mit der „Inneren Energie «“ (nach Blasius „Energie- 
inhalt“) sowie darauf, daß man in der Technik 
unter Strom allgemein eine Größe bezogen auf die 
Zeiteinheit versteht.“ 


Dresden. W. Pauer VDI. 401 
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